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hypothesized,  from  which  false  and  miss-detection 
probability  expressions  are  derived.  Criteria  are  then 
established  for  the  selection  of  sequences  suitable  for  use 
as  unique  words.  Representative  unique  words  of  length  15 
through  20  are  included,  along  with  an  explanation  of  the 
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I.   INTRODUCTION 

The  possible  use  of  time-division  multiple-access  (TDMA) 
in  satellite  communications  has  received  a  considerable 
amount  of  attention  in  recent  years.  In  a  TDMA  system,  many 
stations  use  the  same  satellite  repeater  on  a  time  shared 
basis.  Transmitter  timing  must  be  synchronized  so  that  only 
one  station's  signal  accesses  the  repeater  at  any  given 
time;  and,  further,  so  that  successive  signal  bursts  are 
closely  spaced  in  order  to  make  efficient  use  of  the 
repeater  capacity.  Continuously  changing  ranges  from  ground 
stations  to  the  satellite  (Appendix  B)  make  transmitter 
synchronization  a  non-trivial  problem.  Once  the 
synchronization  problem  is  solved,  TDMA  offers  several 
advantages  over  current  f reguency-division  multiple-access 
procedures  [1/2]o  Schmidt  [3]  outlined  the  proposed 
application  cf  TDMA  to  the  INTELSAI-4  series  commercial 
communication  satellites. 

A.   TDMA  SYSTEM  DESCRIPTION 

As  each  station  must  transmit  in  timed  bursts, 
transmissions  are  necessarily  digital  and  provisions  must  be 
made  for  the  sampling  (and  multiplexing)  of  the  transmitted 
information  channel (s).  This  process  must  be  timed  so  that 
it  is  in  consonance  with  the  burst  repetition  freguency 
(BRF)  and  transmitted  data  rate.  When  analog-to-digital 
conversion  is  involved,  the  analog  sampling  rate  must  be 
egual  to,  or  a  multiple  of,  the  BSF;  a  sample  storage 
capability  is  required  when  the  sampling  rate  is  a  multiple 
of  the  BRF.  When  information  to  be  transmitted  is  in  the 
form  of  a  continuous  bit  stream,  a  buffer  is  reguired  to 
store  a  given  portion  of  the  bit  stream  for  transmission  at 
the  (higher)  data  transmission  rate.  In  either  case,  the 
number  of  information  channels  that  can  be  transmitted  in  a 
burst   of   given  length  is  a  direct  function  of  the  ratio  of 


transmission  rate  to  information  rate.  Time-division 
multiple-access,  then,  is  in  effect  a  higher  level  form  of 
time-division  multiplexing.  TDMA  is  a  time  sharing  at  RF  , 
while  TDM  is  a  time  sharing  at  baseband,  before  carrier 
modulation.  The  process  of  several  transmitted  signals  time 
sharing  a  satellite  repeater  is  conceptually  similar  tc  that 
of  several  information  channels  time  sharing  the  same  RF 
carrier. 

Given  that  transmitted  signals  are  neatly  interleaved  in 
their  arrival  at  the  satellite,  they  will  be  retransmitted 
in  the  same  fashion  and  received  on  the  ground  without  any 
interference  overlap.  A  receiver  must  first  recognize  when 
a  burst  is  of  interest  and  then,  in  the  case  of  multiplexed 
transmissions,  select  particular  information  channels  within 
the  burst  and  route  them  to  appropriate  end  terminals.  In 
order  to  do  so,  the  receiver  must  be  provided  with  certain 
control  information. 

1 •   Receiver  Funct io^a^  Requirements 

Successive  received  bursts  have  been  originated  by 
different,  in  general  widely  separated,  stations  whose 
carrier  frequencies  are  not  coherent  in  frequency  and  phase. 
Therefore,  prior  to  the  demodulation  of  each  burst,  the 
receiver  local  oscillator  must  be  locked  in  frequency  and 
phase  to  that  of  the  incoming  carrier  using  phase-locked 
loop  techniques.  (Use  of  differential  phase-shift  keying 
removes  the  carrier  recovery  requirement  at  the  expense  of 
error-rate  performance  £4].)  Similarly,  once  carrier 
recovery  has  been  accomplished,  the  receiver  bit  timing 
oscillator  (clock)  must  be  synchronized  to  the  incoming  bit 
stream.  Bitzer  (j  5  ]  proposed  a  system  design  for 
establishing  coherent  bit  timing  among  all  stations  and 
Suguri,  Doi  and  Metzger  [4]  reported  on  the  successful  test 
of  this  type  of  system.  In  such  a  system,  receiver  and 
transmitter  clocks  are  continuously  adjusted  to  maintain 
coherency  within  acceptable  limits,  obviating  the  need  for 
gross  adjusteents  to  the  receiver  clock,  preparatory  to  burst 
demodulation. 


With  the  local  oscillator  and  clock  synchronized, 
demodulation  of  the  incoming  signal  occurs.  Embedded 
somewhere  in  the  beginning  of  the  burst  must  be  information 
which  the  receiver  can  extract  to  obtain  transmitting 
station  identity,  word  (channel)  synchronization  and  channel 
status.  Demodulation/demultiplexing  of  the  information 
channels  commences  thereafter. 
2.   l£JJ2S  §•}£  pur st  format 

Satellite  repeater  access  time  is  divided  into 
frames,  a  frame  bein/j  the  length  of  time  allotted  for  the 
transmission  of  one  burst  from  every  station.  The  frame 
length  is  typically  125  usee,  or  a  multiple  thereof,  to 
conform  to  the  Nyguist  sampling  rate  for  voice  (8KHZ)  [3,6]. 
A  station  is  assigned  a  time  slot  within  the  frame  for  the 
transmission  of  its  burst,  the  length  of  which  depends  on 
the  anticipated  channel  requirements.  A  guard  time  of 
typically  200  nsec  or  less  separates  the  time  slcts  to 
preclude  burst  overlap. 

Each  burst  commences  with  a  preamble,  followed  by 
the  information  channels.  A  typical  format  is  depicted  in 
Fig.  1-1.  A  fixed  pattern  of  dummy  bits  is  inserted  at  the 
beginning  for  the  purpose  of  locking  up  the  local  oscillator 
and  clock.  A  constant  string  such  as  0000«®«  for  carrier 
recovery,  followed  by  an  alternating  pattern  such  as  1010««« 
for  bit  timing  recovery  is  typical.  Following  the  recovery 
bits  is  a  fixed  sequence  of  bits  unique  to  the  transmitting 
station,  the  "unique  word."  More  will  be  said  about  the 
unigue  word  below.  .The  remainder  of  the  preamble  contains 
control  signals  and  orderwire  transmissions  necessary  for 
assignment  and  control  of  the  information  channels. 
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Figure  1-1 
Frame  ana  Burst  Format 


TDMA  efficiency  is  a  measure  of  the  total  time 
allotted  to  the  transmission  of  actual  information  versus 
the  total  available  repeater  time.  The  preamble  is  wasted 
transmission  time  in  this  sense  and  considerable  effort  has 
been  expended  to  reduce  its  length  [4,5,8,9], 
3 •   Ujtklcj ue  W or d_  Functions 

lo  the  extent  that  it  is  recognizable  tc  the 
receiver,  the  unique  word  can  be  used  to  identify  the 
originating  station.  It  may  also  be  used  as  a  timing 
reference  to  establish  word  synchronization,  resetting  a  bit 
counter  which  properly  routes  the  control  and  orderwire 
portions  cf  the  preamble  and  indicates  the  beginning  cf  the 
first  information  channel. 

It  has  been  assumed  up  to  this  point  that  all 
transmitting  stations  are  synchronized  so  that  no  burst 
overlap  occurs.  Nothing  has  been  said  about  how 
synchronization  is  acquired  or  maintained.  The  usual  method 
of   doing   sc   is   to   use   the  unique  words  as  timing  marks 
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within  the  frame.  One  station  is  designated  as  the  naster 
station  and  its  unique  word  used  as  the  fundamental  timing 
reference  for  the  frame.  All  other  stations  compare  the 
time  difference  between  their  own  unigue  word  and  the 
reference,  as  received  from  the  satellite,  and  adjust 
transmitter  timing  as  necessary,  one  bit  per  frame,  to 
maintain  prober  separation.  This  problem  of  transaitter 
synchronization  is  of  interest  herein  only  in  that  the 
unigue  word  plays  a  vital  role,  and  will  not  be  further 
pursued.  references  3,7,9  and  10  offer  a  detailed 
description  cf  the  various  factors  which  must  be  considered 
and  provide  insight  as  to  how  a  synchronization  package 
might  be  instrumented. 

B.   DETECTION  AND  SELECTION  OF  UNIQUE  WORDS 

A  TDMA  system  of  the  type  outlined  above  can  be 
operated  successfully  only  if  the  receiver  is  able  to  detect 
the  receipt  of  a  unigue  word  and  differentiate  one  from  the 
other.  Some  type  of  correlation  detector  is  inferred. 
Proposed  detector  design  features  and  operation  are 
presented  in  Section  II.  Then,  in  Section  III,  criteria  are 
established  for  the  selection  of  unigue  words  which  yield 
acceptable  miss-detection  and  false  detection  probabilities. 
Interest  in  the  subject  was  motivated  primarily  by  the 
discussions  of  Schrempp  and  Sekimoto  [11].  Much  of  the 
detector  design  philosophy  is  borrowed  from  their  work,  and 
their  approach  to  unigue  word  selection  has  been 
instrumental  in  obtaining  an  appreciation  for  the 
fundamental  reguireaents  which  must  be  considered.  The 
selection  criteria  set  forth  in  Section  III  and  the 
subseguent  computer  search  results  represent  the  major 
products  cf  this  research  effort. 
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II.   DETECTION  OF  .UNIQUE  WORDS 

A  unique  word  is  a  fixed  binary  sequence  of  length  N. 
Given  a  knowledge  of  the  elements  of  the  sequence,  detection 
can  te  accomplished  .by  comparing  received  bits,  N  at  a  time, 
with  a  replica  of  the  sequence  stored  in  the  receiver;  i.e., 
by  using  correlation  techniques.  The  probability  of  missing 
the  unique  wcrd  (miss-detection)  and  the  probability  of 
receiving  the  sequence  at  random  (false  detection)  are 
fundamental  considerations  in  the  evaluation  of  detector 
effectiveness. 

A.   BASIC  DETECTOR  DESIGN 

A  digital  correlation  detector  of  the  type  depicted  in 
Fig.  2-1  is  postulated.  The  element  values  of  the  unique 
word  are  set  semipermanently  in  a  register  and  the 
demodulated  bit  stream  is  clocked  through  the  shift 
register.  The  corresponding  elements  of  the  two  registers 
are  added,  mcdulo-2,  once  each  bit  period  and  the  results 
are  summed  algebraically  for  input  to  the  threshold  network. 

The  mcdulo-2  operation  is  defined  by, 

0  ffi  0  =  0 

0  9  1  =  1 

1  0  0  =  1 
1  ©  1  =  0 

And  the  input  to  the  threshold  network  is, 


b  =  2   (s  e  u  )  (2.1) 

i=1   i   i 


14 


s 


|S  I        .     .     .  SHIFT    REGISTER  ...     I        S 

1       i  2       I  , >  \  N 

.-. ..,->.,,    .  . ,. .  .  ...    , .  ,  .  .,  .  _^_ ,  .1  __„  ,. .  .  _j 


I  i 

i     1       I     1    2  1     |     N 

U 

4; 

["threshold  1 

{ liii J 

I 

DETECTION 
PULSE 


Figure  2-1 
Digital  Correlation  Detector 


R  therefore  equals  the  number  of  disagreements  between  S 
and  u* .  The  threshold  network  is  designed  so  that  a 
detection  pulse  is  produced  whenever  there  are  h  or  fewer 
disagreements;  that  is,  whenever  R<h.  The  setting  of  a 
threshold  permits  the  detection  of  unique  words  which 
contain  up  to  h  transmission  errors  and  thereby  reduces  the 
miss-detection  probability.  The  probability  of  false 
detection  increases  in  doing  so,  however,  due  to  the 
increased  opportunities  for  a  random  match-up. 

It  should  be  noted  that  R,  as  mathmatically  defined,  is 
in  fact  the  inverse  correlation  function  of  S  and  U.  When 
the  elements  of  S  and  U  all  match,  R=0;  and  when  there  is 
complete  disagreement,  R=N.  Alternative  detector  designs 
can  be  envisioned  which  would  produce  the  true  correlation 
function  at  R;  e.g.,  construct  adders  that  produce  the 
complement  of  the  modulo-2  sum  and  design  the  threshold 
network  so  that  a  pulse  is  produced  whenever  R>N-h.  Of 
course,  in  the  end,  the  voltage  levels  associated  with  '0' 
and   '1'   are   what   determine   the   actual   "value"   of   R. 
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Whichever  of  the  various  alternatives  are  chosen,  the  end 
result  is  the  same:  a  detection  pulse  is  produced  whenever 
there  are  h  or  fewer  disagreements  between  S  and  U . 

The  capability  of  receiving  more  than  one  transmitting 
station  is  provided  by  a  parallel  bank  of  detectors,  each 
set  to  the  unique  wprd  of  one  of  the  transmitting  stations. 
The  received  bit  stream  is  passed  simultaneously  through 
each  detector  and  the  resulting  detection  pulses  are  used  to 
control  receiver  functioning,  as  described  in  Section  I. A. 3. 

B.   STEADi-STATE  APERTURE  GATING 

The  probability  of  false  detection  is  reduced 
significantly  if  the  detector  can  be  gated  so  that  it  begins 
looking  for  a  unigue  word  just  prior  to  its  receipt,  and  is 
otherwise  inactive.  The  detection  pulse  itself  may  be  used 
to  generate  such  a  gating  signal  under  steady-state 
conditions;  that  is,  once  the  unigue  word  has  been  initially 
acguired . 

A  unigue  word  it*  transmitted  every  I  seconds,  once  each 
frame,  except  when  a  correction  is  reguired  to  maintain 
proper  position  within  the  frame.  If  corrections  are  made 
at  the  rate  of  one  bit  per  frame  [9,10],  the  elapsed  time 
between  unigue  word  transmissions  is  bounded  by  T±b,  where  b 
is  the  period  of  one  bit.  The  time  separation  of  successive 
received  unigue  words  is  further  affected  by  relative  motion 
between  the  satellite  and  the  respective  transmitting  and 
receiving  stations.  Let  d,  in  bits  per  second,  be  the 
maximum  expected  doppler  induced  drift  of  successive  unigue 
word  detection  pulses.  The  caximum  doppler  shift  in  one 
frame  period  is  then  ±Tbd  and  the  elapsed  time  between 
successive  receipts  of  a  unigue  word  is  bounded  by, 

T   =  T  ±  (b  +  Tbd) 
r 

T   =  T  ±  b(1  +  Td)  (2.2) 

r 
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This  relationship  can  be  expressed  in  terms  of  bit 
periods  by  dividing  through  by  b.  Doing  so  establishes  a 
limit  on  the  number  of  bit  periods  separating  successive 
receipts  of  a  unique  word. 

N   =  T  /b  =  T/b  ±  (1  +  Td) 
r    r 

N   =  n  ±  (1  +  Td) 

r 

where  n  is  the  number  of  bits  in  one  frame. 

Defining  D  as  tbe  maximum  displacement  {in  bits)  cf  the 
unique  word  from  its  "periodic"  position, 

D  =  1  +  Td  (2.3) 

N   =  n  ±  D  (2.4) 

r 

Upon  receipt  of  a  unigue  word,  the  detector  is  turned 
off,  gated  back  on  T-bD  seconds  later,  gated  off  by  the  next 
unique  uoid,  etc.  False  detections  can  therefore  result 
only  from  a  random  match-up  within  the  aperture  immediately 
preceding  the  unigue  word. 

The  guestion  remains  as  to  what  to  do   when   the   unigue 

word   is   missed   due   to   errors  occurring  in  transmission. 

Once  the  detector  is   gated   on,   is   it   left   on   until  a 

detection   pulse   is  received  (in  which  case  a  missed  unigue 

word  causes  an   immediate   transistion   to   the   acquisition 

mode) ;   or,  is  the  detector  gated  off  at  time   T+bD  and  then 

1 
gated  back  on  at  time   2T-bD,   one   frame   later?     In   the 

latter  case,   how  many  successive  misses  are  required  before 

cancelling  the  aperture   and   reverting   to   an   acquisition 

mode?    The   answer   depends  for  the  most  part  on  particular 


1 

It   is  assumed  that  the  aperture  width  remains  constant 
and  would  not  be  expanded  as  successive  misses  occur. 
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system  specifications.  When  a  unique  word  is  missed,  all 
information  contained  in  the  burst  is  lost.  Having  missed 
it  once,  that  is  having  received  the  unigue  word  with  more 
than  h  errors,  there  is  a  given  probability  that  it  will 
contain  more  than  h  errors  the  next  time.  Calculation  of 
the  probability  of  a  unigue  word  containing  more  than  h 
errors  two  or  more  times  in  succession  is  not  enough 
however.  The  doppler  induced  burst  drift  rate  must  also  be 
considered. 

If  the  drift  rate  is  such  that  the  unigue  word  is 
received  outside  tiie  aperture,  it  will  be  missed,  even  if 
error  free.  On  the  pther  hand,  reversion  to  an  acguisition 
mode  may  involve  a  considerable  penalty  in  lost  information 
due  to  a  long  average  acguisition  time.  (See  Section  II. C 
below.)  Thus,  a  desire  to  avoid  the  acguisition  mode  must 
be  balanced  against  the  diminishing  probability  (as 
successive  misses  occur)  of  receiving  the  unigue  word  in  the 
aperture  mode.  As  stated  earlier,  this  is  a  systems 
dependent  decision,  and  one  that  can  best  be  reached 
empirically  in  an  operational  environment.  For  present 
purposes,  the  rather  safe  assumption  will  be  made  that  at 
least  two  successive  misses  are  required  before  a  transition 
from  the  aperture  to  acguisition  mode  is  effected.  A 
detection  pulse  turns  off  the  detector  and  sets  a  counter 
which  gates  the  detector  on  T-bD  seconds  later.  If  no 
detection  pulse  occurs,  the  counter  turns  off  the  detector 
at  time  T+bD.   The  aperture  is  of  width  2D-H  (bits) . 

The  maximum  burst  drift  which  may  be  expected  when  using 
a  geosynchronous  satellite  is  computed  in  Appendix  B. 
Eguation  B.5  can  be  substituted  into  eguation  2.3  to  obtain 
values  for  the  maximum  displacement  of  a  unigue  word  over 
one  frame  period.  Table  2-1  contains  representative 
calculated  results  for  the  shipboard  platform  discussed  in 
Appendix  B  as  well  as  for  a  platform  with  a  relative 
velocity  on  the  order  of  1000  knots. 
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Table  2-1 
Bit  Rate  Doppler(d)  and  Maximum   Displacement (D) 
of  a  Unique  word  Over  One  Frame  Period  for  Given 
Values  or  Station/Satellite  Relative   Motion  (v)  , 
Bit  Frequency (f)  and  Frame  Period (T) 


C.   ACQUISITION  GATING 

The  steady-state  operation  of  a  unique  word  detector  was 
described  in  the  preceding  section.  Here  the  process  of 
acquisition  is  considered.  When  a  receiver  is  activated, 
each  detector  must  acquire  the  associated  unique  word  within 
the  TDMA  frame.  Likewise,  a  detector  must  be  switched  to" 
the  acquisition  mode  whenever  steady-state  synchroni2ation 
is  lest.  Acquisition,  or  re-acquisition,  may  be  effected  by 
simply  passinq  the  received  bit  stream  through  the  detector 
until  a  detection  pulse  is  received.  Then,  as  in  the 
steady-state  mode,  the  detection  pulse  switches  the  detector 
off  and  starts  an  aperture  qate  counter  which  in  turn 
switches  the  detector  back  on  in  the  next  frame.  If  a 
second  detection  pulse  is  received  within  the  aperture,  the 
detector  is  transferred  to  the  steady-state  mode. 
Otherwise,  the  acquisition  process  continues. 

When  in  the  acquisition  mode,  the  bursts  of  all  active 
stations  are  being  passed  continuously  through  the  detector 
and  the  probability  of  false  detection   is   thus   relatively 
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high.  Should  a  false  detection  occur,  a  false  aperture  is 
set  and  the  next  actual  unique  word  detection  pulse  will  be 
missed.  Even  when  the  detection  is  valid,  there  is  a 
non-zero  probability  that  the  unigue  word  will  be  missed  in 
the  next  frame.  la  either  case,  acquisition,  and  transfer 
to  the  steady-state  mode,  is  delayed  for  at  least  two  frame 
periods.  The  actual  acquisition  time  may  thus  span  several 
frame  periods. 

It  is  possible  to  construct  a  detection  pulse  storage 
network  which  would  look  at  detection  pulse  locations  over 
two  or  more  frame  periods  and  use  majority  logic  to  set  an 
aperture  corresponding  to  the  true  detection  pulse. 
Although  acquisition  time  would  be  improved,  implementation 
could  prove  relatively  costly  because  each  detector  would 
reguire  its  own  network.  An  alternative  is  proposed,  which 
makes  use  of  the  same  aperture  network  employed  in  the 
steady-state  mode  to  gate  the  detector  and  reduce  the 
probability  of  false  detection. 

1  •  lE§:®§   Sefg^eijce  Gating 

Each  burst  has  an  assigned  position  within  the  TDMA 
frame.  It  is  known  that  a  particular  burst  will  be  received 
a  given  number  of  bits  after  the  receipt  of  the  reference 
unigue  work,  plus  or  minus  the  guard  band  cushion.  "  The 
guard  band  is  set  so  that  it  prevents  burst  overlap;  that 
is,  it  represents  the  maximum  expected  displacement  of  a 
burst  (or  unigue  word)  from  its  nominal  position.  It  is 
therefore  proposed  that,  during  acquisition,  the 
aperture/counter  network  of  a  detector  be  keyed  by  the 
receipt  of  the  reference  unique  word.  The  aperture  is  set 
so  that  the  detector  is  gated  on  at  the  assigned  position  of 
the  corresponding  unique  word,  minus  the  number  of  bits  in 
the  guard  band,  B.  If  the  unigue  word  is  missed,  the 
detector  is  switched  off  B  bits  after  the  nominal  unique 
word  position.  (The  aperture  width  is  2B+1.)  If  the  unique 
word  is  detected,  the  aperture  is  set  for  the  next  frame. 
Two  consecutive  detections  key  a  transfer  to  the 
steady-state  mode. 
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In  a  demand  access  system,  where  at  least  some  of  the 
stations  are  not  assigned  a  permanent  time  slot,  use  of 
frame  referecce  gating  is  feasible  only  if  all  stations  can 
be  provided  with  time  slot  activation  and  deactivation 
information,  or  can  otherwise  predict  where  a  given 
station's  burst  may  appear  within  the  frame.  Whether  cr  not 
it  is  possible  to  provide  advance  activation  notification  to 
all  stations  depends  on  how  a  part  time  station  gains  access 
to  a  time  slot.  If  access  is  centrally  controlled 
(automatically) ,  then  assignment  signals  from  the 
controlling  station  can  be  copied  by  all  stations  and  used 
to  set  the  acquisition  gate.  If  a  station  gains  access  by 
seizing  a  vacant  time  slot,  a  receiving  station  must  have 
other  means  of  predicting  where  the  burst  will  occur.  The 
results  of  the  computer  search  for  unique  words,  reported  in 
Section  IV,  suggest  a  procedure  for  meeting  this 
requirement. 

For  a  given  unique  word  length,  several  different  sets 
are  available  which  .have  good  cross-correlation  properties. 
Unigue  words  from  any  one  set  could  be  assigned  to  certain 
stations  which  collectively  have  exclusive  use  of  a  block  of 
consecutive  time  slots  within  the  frame.  As  many  such 
blocks  as  necessary  could  be  made  available  to  different 
groups  of  stations*.  The  number  of  time  slots  assigned  to 
each  group  wculd  be  based  on  collective  useage  requirements, 
and  could  be  varied  on  a  scheduled  or  demand  basis.  (Demand 
access  implies  that  the  number  of  time  slots  would  be  less 
than  the  number  of  stations.)  Those  receiver  detectors 
corresponding  to  the  demand  access  transmitters  could  then, 
when  in  the  acquisition  mode,  be  gated  on  at  the  beginning 
of  each  burst  withi-c  the  demand  access  sub-frame.  The 
probability  of  false  detection  will  be  somewhat  higher  than 
for  the  single  gate  procedure  proposed  above;  but,  with  low 
cross-correlation  between  unique  words  of  a  set,  will  be 
much  lower  than  for  ."wide  open"  acquisition. 
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Ill-   BNIQPE  WORD  SELECTION  CRITERIA 

Having  outlined  the  basic  detector  design  concepts,  a 
basis  has  been  formed  for  the  selection  of  sets  of  unique 
words.  A  detailed  analysis  of  the  attendant  detection 
probabilities  is  presented  in  Section  C.  Simply  stated,  the 
goal  is  to  find  unique  words  which  yield  low  miss-detection 
and  false  detection  probabilities.  Criteria  will  be 
established  in  this  section  for  doing  so. 

A.   LENGTH 

The  fewer  bits  in  a  unique  word,  the  lower  the 
miss-detection  probability.  Conversely,  the  longer  the 
unigue  word,  the  lower  the  false  detection  probability.  The 
use  of  a  detector  threshold  and  aperture  nullifies  to  some 
extent  the  importance  of  these  opposing  considerations,  but 
in  a  gross  sense  the.y  can  never  be  eliminated  as  factors 
influencing  the  ultimate  choice  of  unique  words.  In  fact, 
length  and  the  threshold  setting  are  the  only  parameters 
available  for  altering  the  raiss-detection  probability. 

From  the  standpoint  of  efficiency,  it  is  desirable  to 
keep  unigue  words  as  short  as  possible.  The  longer  each 
station's  unigue  word,  the  less  time  is  available  for  the 
transmission  of  information. 

E.   AUTOCORRELATION  CRITERIA 

Operation  of  the  detector  in  the  steady-estate  mcde  is 
described  in  Section  II. B.  Loss  of  channel  information  will 
result  if  (1)  the  unique  word  is  missed,  or  (2)  a  false 
detection  occurs  in  the  aperture  prior  to  the  unique  word. 
Either  may  lead  eventually  to  a  drop-sync  condition.  As  (1) 
and  (2)  are  equally  damaging,  it  would  seem  appropriate  to 
select  unique  words  in  such  a  way  that  their  combined  effect 
is  minimized.  Under  the  conditions  described  in  Section 
II.fi,   it   is   largely   a   matter   of   unquantif iable  chance 
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whether  a  miss-detection  or  a  false  detection  (within  the 
aperture)  is  the  least  desirable  occurrence.  Accordingly, 
the  following  condition  was  established:  if  the  received 
unigue  word  contains  h  or  fewer  errors,  the  false  detection 
probability  must  be  zero  within  the  aperture.  A  false 
detection  pulse  will  not  t>e  produced  unless  the  unigue  word 
is  going  to  be  missed  anyway. 

As  it  is  possible  for  a  unigue  word  detection  pulse  to 
occur  on  the  last  cycle  within  the  aperture,  and  as  the 
aperture  is  of  width  2D*1,  the  foregoing  condition  can  be 
met  only  if  the  probability  of  false  detection  is  zero  for 
the  2D  cycles  preceding  the  unigue  word.  This  condition  can 
in  turn  be  satisfied  only  if  the  autocorrelation  function  of 
the  unigue  word  is  such  that, 

fi;(t)  <  N  -  2h  -  t  ;   1<t<2D       (3.1) 

where  t  is  an  integral  time  variable,  the  basic  unit  of 
which  is  cne  bit  period,  and,  . 

N-t 
fi(t)  =.2  U(i)©U(i  +  t) 

i=1 

This  condition  is  presented  graphically  in  Fig.  3-1. 
R  (t)  is  egual  to  the  number  of  agreements  between  the 
(error-free)  N-t  displaced  unigue  word  bits  and  the  detector 
setting.  Up  to  h  errors  are  permitted  in  the  unigue  word 
under  the  prescribed  conditions  and  it  is  possible  that 
these  errors  may  lie  such  as  to  produce  additional 
agreements,  giving  a  total  of  R  (t) +h  agreements  frcm  the 
displaced  unigue  word  alone.  There  are  t  non-unigue  word 
bits  in  the  detector.  If  all  of  them  match  with  the 
detector  setting,  a  combined  total  of  R(t)+h  +  t  agreements 
results.  As  N-h  or  more  agreements  are  reguired  tc  produce 
a  false  detection,  the  probability  of  such  occurrence  is 
zero  if  E  (t) +h*t<N-h,  and  the  validity  of  ineguality  3.1  is 
established. 
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In  the  following  section,  a  condition  for  the 
cross-correlation  function  will  be  established.  This  same 
condition  must  be  met  by  the  autocorrelation  function 
outside  the  aperture  in  order  to  meet  specified  limits  on 
the  probability  of  false  detection  during  acquisition.  It 
is  listed  here,  in  advance,  in  order  to  complete  the 
autocorrelation  specifications. 

R  (t)  <  N  -  h  -  t  ;   2D<t<N-h        (3.1a) 

£(t)  =  P  ;   N-h<t<N-1       (3.1b) 

N 
There  are  2   seguences  of  length  N  from  which  to   choose 

unigue   words.     However,   symmetry   considerations   can  be 

invoked  to  reduce  the  number  which  must  actually  be  checked. 

First,  because  the  mpd-2  summing  operation  is  commutative,  a 

seguence  and  its  complement,  obtained  by  inverting  the  value 

of   each   bit,   have   the   same   autocorrelation   function. 

Second,  because  autocorrelation   is   an   even   function,   a 

sequence   and   its   transpose   (or   mirror) ,   obtained   by 

inverting  the  order  of  the  bits  in  the   seguence,   have   the 

same   autocorrelation   function.   Therefore,  once  a  seguence 

is  checked   it   is   -unnecessary   to   check   its  complement, 

transpose   or  complement   of   its  transpose.   The  number  of 

reguired  checks  is  npt  reduced  by  a  factor  of  four,   as   it 

would  appear  at  first  glance,  however,  because  the  transpose 

and  complement   transpose   pair   may   be   identical   tc   the 

original   seguence   and  its  complement.   The  reguired  number 

of  checks  can  be  determined  by  consideration  of   the   effect 

of   the   complement   and  transpose  operations.   Although  the 

results  were   not   used   directly   in   the   computer   search 

routine,  they  do  provide  an  indication  of  search  efficiency, 

and  are  presented  here  as  a  matter  of  peripheral  interest. 

In  order  for  the  transpose  or  complement  transpose  to  be 

identical  to  the  original  seguence,  the  second  [N/2]  bits  of 

the   seguence   must  be  the  transpose  or  complement  transpose 

of  the  first  [N/2]  bits,   where   [N/2]   equals   the   largest 
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integer  less  than  or  equal  to  N/2.  (If  N  is  odd,  the  middle 
bit  of  the  sequence  does  not  change  position  under  the 
transpose   operation,.)   The  first  [N/2]  bits  may  be  arranged 

[N/2] 

in   2       different   ways,   each   arrangement  having   an 

associated   transpose  and  complement  transpose  in  the  second 

[ N/2  ]  N 

half.   Therefore,  2x2      of  the  original  2   sequences  have 

redundant   (identical   to   original   sequence)  transposes  or 

complement  transposes,  and,  for  this   subset,   the   required 

number   of  checks  can  be  reduced  by  a  factor  of  two  only  (by 

pairing   a   sequence   with   its   complement)^     There   are 

N   [N/2]+1 
2  -2         sequences  remaining  to  be  considered.   These  can 

be   grouped   by   fours   (sequence,   complement,    transpose, 

complement   transpose)   to   reduce  the   number   of  reguired 

checks  by  four.   The  total  number  of  sequences  of   length   N 

which   must  be  checked  in  search  of  unique  words  meeting  the 

autocorrelation  condition  is  therefore, 

N     [N/2 1+1 
2-2  £N/2J 

Required  Checks  =  > +  2 

4 

N-2     £  N/2  ]-1 
=  2     +2 

The     computer    program    used    to    accomplish    the 
autocorrelation  sort  is  listed  in  Appendix  D. 

C.   CROSS-COfiRELATIGN  CRITERIA 

The  cross-correlation  between  the  unique  words  used  by 
stations  sharing  the  TDMA  frame  is  a  crucial  factor 
affecting  the  false  detectic  probability  during  (ungated) 
acquisition.  A  cross-correlation  characteristic  that  yields 
a  high  probability  of  false  detection  is  particularly 
unsatisfactory  because  of  the  periodicity  of  its  occurrence. 
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In  such  cases,  the  detection  aperture   may   be   falsely   set 
time  and  time  again  over  several  consecutive  frames. 

The  probability  Expressions  derived  in  Appendix  C  were 
used  as  a  basis  for  the  establishment  of  the  following 
cross-correlation  cpndition:  the  probability  of  false 
detection  must  be  zero  unless  the  received  unique  word 
contains  cne  or  more  errors.  This  condition  also  applies  to 
the  head  and  tail  of  the  autocorrelation  function,  outside 
the  aperture.  The  cross-correlation  condition  is  thus 
similar  to  the  autocorrelation  condition,  differing  only  in 
that  the  unigue  word  must  be  error  free,  vice  h  or  fewer 
errors,  to  guarantee  zero  false  detection  probability. 
Using  fi(t)  again  to  xepresent  the  correlation  function,  th<e 
condition  en  R (t)  is  (analogous  to  ineguality  (3.1)) 
R  (t) <N-h- 1 1 1 .  As  the  correlation  function  is  always 
non-negative,  this  condition  can  be  met  only  when  |t|<N-h. 
For  higher  values  of  jt|  the  condition  must  be  relaxed  to 
R(t)=0.  Summarizing,  sets  of  possible  unique  words,  having 
passed  the  autocorrelation  test,  are  retained  only  if  the 
cross-correlation  between  all  words  of  the  set  satisfy, 

R(t)  <  N  -  h  -  |t|  ;   |t|<N-h       (3.2a) 

R(t)  =  0  ;   N-h<|t|<N-1   (3.2b) 

The  conditions  placed  on  R  (t)  can  be  used  to  further 
refine  the  probability  expressions  of  Appendix  C.  It  will 
be  assumed  that  bits  of  fixed  value  (at  least  N-1  of  them) 
precede  the  unique  word  and  that  bits  of  random  value 
follow.  This  assumption  is  based  on  the  system 
specifications  outlined  in  Section  I,  and,  although  not  an 
absolute  necessity  at  this  point,  is  made  in  order  to  reduce 
the  number  of  situations  which  must  be  considered.  Other 
situations,  when  and  if  occurring,  may  be  handled  in  a 
completely  analogous  fashion.. 

First,  the  probability  of  false  detection  in  the  fixed 
bit  region  preceding  a  unique  word  (a  unique  word  other  than 
the  one  for  which  the  detector  is  set)  will   be   considered. 
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The  point  where  the  unique  word  is  completely  contained  in 
the  detector  defines  t=0.  At  t=-(N-1)  the  first  bit  cf  the 
unigue  word  enters  the  detector.  Eguation  (s)  C.8  apply  for 
-(N-1)<t<0,  with  m=N.  As  in  Appendix  C  "number  of 
agreements11  refers  to  matches  between  the  detector  setting 
and  the  value  of  the  fixed  bits  as  transmitted.  The  number 
of  agreements,  k  in  eguation  C.8,  is  equal  to  R (t)  ,  the 
number  of  matches  within  the  unique  word  bits,  plus  the 
number  of  matches  within  the  remaining  jt|  bits.  If  all  of 
the  latter  match, 

k  =  R  (t)  <-  |t  J 

Therefore,  the  cross-correlation  condition  3.2a 
(R  (t)  <N-h- 1 1  J)  is  equivalent  to  the  equation  C.8  condition, 

k  +  h  <  N 

The  cross-correlation  condition  3.2b  in  itself  does  not 
guarantee  that  k+h<N.  It  is  necessary  to  impose  the 
additional  systems  dependent  condition  that  fewer  than  K-h 
agreements  can  be  permitted  between  the  non-unique  word  bits 
for  -  (N-1)  <t<- (N-h)  .  Doing  so,  the  probability  of  false 
detection  at  any  given  step  is  (from  eguation  C.8e), 

P{A<t)}  = 

N-h-k      h+k 
C(N-k,N-h-k)  p      (1-p)    ;   -(N-1)<t<0    (3.3) 

If  k  takes  on  its  maximum  value  (N-h) ,  the  probability 
of  false  detection  becomes, 


N-1 
P{A(t)}     =  (h+1)p(1-p)    ;   -(N-1)<t<0   (3.4) 
max 


The  probability  of  false  detection  in  the  random  bit 
region  following  the  unique  word  is  also  given  by 
eguation (s)  C.8.   In  this  case  m  equals  the  number  of  unique 
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word  bits  in  the  detector  (there  are  no   other   fixed   bits) 
and  the  number  of  fixed  matches,  k,  equals  fi  (t) . 

m  =  N  -  |t| 

k  =  E(t) 

Therefore,    the   cross-correlation   condition    3.2a 
(R  (t) <N-h- 1 1 J)  is  equivalent  to  the  equation  C.8  condition, 

k  ♦  h  <  m 

And   the  probability  of  false  detection  at  any  given  step  is 
given  by  equation  C.Bc   (with  m  and  k  given  by  above) . 

P{A(t)}  = 

N-m  m-k-h      k+h 

(1/2)     C(m-k, m-k-h)  p      (1-p)    ;   0<t<N-h    (3.5) 

When   k   takes   on   its   maximum   value   (R (t) =N-h-t-1) , 
equation  3.5  becomes. 


t  N-t-1 

P{A(t)}     =  (1/2)  (h+1)p(1-p)      ;   0<t<N-h    (3.6) 
max 


For  N-h<t<N-1,  m=h,h-1,h-*2,®®<»,1  and  k=R(t)=0,  giving 

m  =  N  -  t 

and, 

k  +  h  >  m 

The  protability  of  false   detection   is   (from   equation 
C.8b)  , 

N-m      m   h-^m 
P{A(t)}  =  (1/2)    (1-p)   5    C(N-m,j)         (3.7) 

3  =  0 

t      N-t   h-N  +  t 
=  (1/2)     (1-p)     5      C(t,j);    N-h<t<K-1 

3  =  0 
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Equation  3.7,  unlike  equations  3.3  and  3.5,  is 
independent  of  parameter  variation  [i.e.  independent  of  k, 
fi  (t)  ]  and  thus  "maxi<mum  P{A(t)}"  has  no  meaning. 

Equations  3.3,  3.5  and  3.7  may  be  used  to  calculate  the 
probability  of  false  detection  as  each  unique  word  in  the 
frame  is  passing  tiarough  the  detector.  Equations  3.4,  3.6 
and  3.7  provide  a  worst  case  limit  for  these  probabilities 
(when  unique  words  are  chosen  usinq  the  cross-correlation 
criteria  3.2a  and  3.2b).  These  individual  fixed-bit  event 
probabilities  are  the  ones  needed,  among  others,  in  equation 
C.10  to  obtain  an  upper  bound  for  the  probability  of  one  or 
more  false  detections  (F)  over  the  period  of  one  frame.  Or, 
over  a  narrower  time  span,  they  can  be  used  to  compute  an 
upper  bcund  on  the  probability  of  one  or  more  false 
detections  durinq  the  2N-1  detection  events  surrounding  a 
given  unigue  word;  that  is,  to  compute, 


N-1 
t=-N+1 


F1  provides  a  means  of  obtaining  a  relative  measure  of 
the  goodness  of  different  unique  word  choices  without  having 
to  resort  to  the  more  lengthy  calculations  required  by 
equation  C.10.  Substituting  from  equations  3.4f  3.6  and 
3.7, 

N-1      N-h      t  N-t-1 

F«  <  N[  (h*1)p(1-p)    ]  +  5    (V2)   (h+1)p(1-p) 

t=1 

N-1        t      N-t   h+t-N 
**?„   t  ,  (V2)  (1-p)     ,2n         c(t,j)         (3.8) 
t=N-h-1  3=0 

The  computer  proqram  used  to  check  cross-correlaticn  is 
listed  in  Appendix  D.  Each  sequence  which  passes  the 
autocorrelation  test  is  cross-correlated  with  every  other 
such  sequence,  yielding  pairs  of  acceptable  unique  words. 
In  crder  to  form  maximal  length  sets  of  acceptable  unique 
words,   it   is   necessary   to  perform  an  additional  matching 
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routine.  For  example,  if  UW1  cross-correlates  well 
(matches)  with  UW2*  UW3  and  UW4,  then  a  set  of  three  is 
obtained  if  UW2  also  matches  with  UW3  or  UWU,  or  UW3  matches 
with  UWU;  a  set  of  four  is  obtained  if  UW2  matches  with  both 
UH3  and  UW4  anj  UW3  matches  with  0W4.  Various  matching 
algorithms  were  programmed,  none  of  which  performed  in  a 
completely  acceptable  fashion.  At  N=17,  the  number  of 
checks  which  must  be  performed  becomes  so  great  that,  for 
the  algorithms  used,  computation  time  runs  into  days.  In 
order  to  ottain  results  for  N=17  and  higher,  various 
expedients  were  employed,  as  described  in  the  following 
section,  to  reduce  the  required  computer  time. 
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IV.   PRESENTATION  OF  RESULTS 

The  SES-9300  digital  computer  was  used  to  perform  a 
correlation  sorting  routine  patterned  after  the  criteria 
established  in  the  preceding  section.  At  the  outset,  it  was 
realized  that  a  significant  amount  of  machine  time  would  be 
reguired  to  obtain  results;  but,  even  so,  the  magnitude  of 
the  profclem  was  not  fully  appreciated  until  program 
execution  was  initiated.  The  SDS-9300  was  chosen  over  the 
IBM-360  because  of  the  local  availability  of  multi-hour 
blocks  cf  computer  t,ime,  realizing  that  its  limited  on-line 
memory  (32K)  might  prove  restrictive. 

The  choice  of  detector  parameters  and  the  essential 
findings  of  the  computer  search  effort  are  covered  in  the 
following  paragraphs..  Selected  sets  of  unigue  words  are 
listed  in  Appendix  A. 

A.   CHOICE  OP  DETECTOR  APERTURE  AND  THRESHOLD 

The  ffotivation  for  use  of  a  steady-state  aperture 
technigue  is  provided  in  Section  II. B.  As  the  aperture 
width  must  be  at  least  2D+1,  where  D  is  the  maximum 
displacement  of  a  unigue  word  in  one  frame  period,  Table  2-1 
leads  one  to  a  choice  of  a  four  or  five  bit  width.  An  odd 
width  is  preferred  ioa  order  that  the  aperture  counter  can  be 
run  directly  off  the  receiver  clock;  i.e.,  so  that  detector 
activaticn/deactivation  occurs  at  the  beginning  rather  than 
in  the  middle  of  a  bit  cycle.  An  aperture  width  of  five 
bits  is  thus  the  mimimum  acceptable,  and  yet  at  the  same 
time  provides  sufficient  margin  to  cover  unigue  word  drift 
over  more  than  one  frame  in  the  event  of  a  miss-detection. 

The  choice  of  a  detector  threshold  (h)  is  determined  by 
the  setting  cf  an  acceptable  limit  for  the  raiss-'detection 
probability  (Q) ,  which  in  turn  is  a  function  of  the  bit 
error  rate  (p) .   For  this  purpose,  p  was  assigned  a  value  of 
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-4 

10     [ 1 1 3  as  being  the  worst-case  bit  error  rate  that  could 

be  tolerated  on  a  satellite  communication  channel. 
Eguations   C.1,   C.2   and  C.3  were  used  to  compute  Q  and  the 

mean  time  between  miss-detections  (H)  for  various  values   of 

N   and   h.    The  res-alts  are  listed  in  Table  4-1.   The  times 

given  for  T!  assume  a  frame  period  of  125  usee. 

Acceptable   values   for  "H  occur  for  h=2  and  higher.   The 

fact  that  the  probability  of  false  detection  increases  with 
h  dictates  the  selection  of  the  minimum  acceptable 
threshold;  viz.  h=2. 


N 

h 

2   -3 

15 

0 

1.50x10 

-6 

1 

1.05x10 

-10 

2 

4.55x10 

-13 

.3 

1.36x10 

-3 

17 

0 

1.70x10 

-6 

1 

1.36x10 

-10 

2 

6.79x10 

-13 

3 

2.38x10 

-3 

20 

0 

2.00x10 

-6 

,1 

1.90x10 

-9 

2 

1.  14x10 

-13 

3 

4.84x10 

1 

2  min 
76  hrs 

29  yrs 
74  msec 
92  sec 
51  hrs 
17  yrs 
62  msec 
66  sec 

30  hrs 
8  yrs 


Table  4-1 
Miss-Detection  Probabilities 
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B.   UNIQUE  WOBD  SYMMETRIES 

Certain  characteristics  of  the  sequences  which  pass  the 
autocorrelation  test  were  used  to  simplify  the 
cross-correlation  sort.  These  same  characteristics  also 
provide  a  degree  of  flexibility  in  applying  the  results  to 
particular  operating  systems. 

A  seguence  can  pass  the  autocorrelation  test  only  if  its 
autocorrelation  is  .zero  for  the  first  and  last  three  values 
at  the  head  and  tail  of  the  distribution  (with  h=2)  . 
Therefore,  the  first  three  bits  of  the  sequence  must  be  all 
zeros  or  all  ones,  and  the  last  three  bits  must  all  be  of 
opposite  polarity.  The  need  to  check  less  than  one^-half  of 
the  total  possible  seguences  was  explained  in  Section  III.B. 
As  the  autocorrelation  sort  looks  at  seguences  in  ascending 
order  of  numerical  value,  seguences  passing  the  test 
necessarily  begin  with  three  zeros.  (The  search  is 
exhausted  prior  to  reaching  seguences  which  have  a  one  as 
the  first  bit  on  the  left.)  Hhen  a  potential  unique  word 
(UW)  seguence  passes  the  autocorrelation  test  it  is  known 
that  its  complement  (CUW)  ,  transpose  (TUfc)  and  complement 
transpose  (C1UH)  will  also  pass  the  test.  The  latter  were 
computed  and  stored  along  with  the  UW.  For  example,  the 
following  seguences  pf  length  14  pass  the  autocorrelation 
test: 

U«  CUW  TUH  CTU^ 

00000101100111  11111010011000  11100110100000  00011001011111 
00001011001111  11110.100110000  11110011010000  000011001C1111 

The  cross-correlation  condition  is  identical  to  that  of 
autocorrelation  on  the  head  and  tail  of  the  distribution. 
That  is,  the  cross-correlation  between  seguences  must  be 
egual  to  zero  for  the  first  and  last  three  values.  It  is 
obvious  from  the  example  above  that  crossi-correlations 
between  the  Uw/CTUW  sequences  and  the  CUH/TUK  sequences  do 
not  meet  this  condition.   This  is  true  in  general.   Further, 
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because  of  the  symmetry  inherent  in  the  complementing 
operation,  (1)  if  UW1  cross-correlates  well  (matches)  with 
CTUW2,  then  CUW1  matches  with  TUW2  and  (2)  if  UW1  matches 
with  UW2,  then  CTUW1  matches  with  CTUW2,  etc.  Therefore,  in 
checking  cross-correlation  it  is  only  necessary  to  check 
each  UW  against  otier  OW's  and  CTUW's.  Each  set  of  unigue 
words  obtained  by  this  process  can  be  converted,  by 
inspection,  to  a  second  set  by  taking  the  cottplement  of  each 
unigue  word.  This  provides  a  dual  flexibility  in  applying  a 
set  to  a  particular  operating  system  with  its  own  peculiar 
arrangement  of  fixed  value  bits  preceding  the  unigue  word. 

C.   UNIQUE  WCPD  SETS 

As  explained  above,  the  seguences  from  the  classes 
designated  CUW  and  TUW  need  not  be,  and  were  not,  considered 
in  the  cress-correlation  search  for  sets  of  unigue  words. 
Each  set  extracted  from  the  UW/CTUW  seguences  has  a 
companion  within  the  CUW/TUW  seguences,  obtainable  by 
inspection.  In  ordei.  to  accurately  reflect  the 
computational  reguirements,  the  reported  results  which 
follow  thus  are  drawn  only  from  the  former.  That  is,  such 
guantities  as  "number  of  seguences  passing  the 
autocorrelation  test."  and  "number  of  sets  of  length  X"  refer 
to  actual  computer  results.  The  numbers  can  be  doubled  by 
including  the  complemented  seguences.  Computer  program 
descriptions  are  included  in  appendix  D. 

Two  results  of  primary  interest  were  obtained  for  each 
value  of  N.  First,  the  unigue  word  which  cross-correlated 
well  with  the  maximum  number  of  seguences  was  recorded, 
along  with  the  corresponding  matching  seguences.  These  are 
referred  to  as  "gated  sets."  These  seguences  may  be  used  as 
unigue  words  in  a  gated  acguisition  system  where  it  is 
necessary  that  only  the  reference  unigue  word 
cross-correlate  well  with  the  other  unique  words.  That  is, 
as  all  detectors  except  the  one  set  for  the  reference  are 
gated,  only  the  latter   must   contend   with   non-associated 
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unique  words  during  acquisition,  and  is  thus  the  only  one 
for  which  cross-correlation  properties  are  of  concern. 
Second,  sets  of  unique  words  were  obtained  which  have  the 
property  that  every  word  in  the  set  cross-correlates  well 
with  every  other  .word  in  the  set.  Sets  of  this  type  are 
necessary  for  ungated  acquisition,  and  are  referred  to 
hereinafter  simply  as  "sets." 

With  h=2  and  a  five  bit  aperture  width,  no  sequences 
passed  the  autocorrelation  test  for  N=13  and  below.  At 
N=14,  fouL  sequences  passed  the  autocorrelation  test,  but 
none  of  them  could  be  paired  in  cross-correlation.  At  N=15, 
36  sequences  passed  the  autocorrelation  test.  A  qated  set 
of  length  13  and  a  total  of  16  sets  of  length  three  were 
obtained.   The  total  computer  search  time  was  25  sees. 

At  N=16,  150  sequences  passed  the  autocorrelation  test, 
and  the  largest  gated  set  was  of  length  66.  A  total  cf  56 
sets  of  length  six  were  found.  The  computer  search  time  was 
3  min,  50  sees. 

The  difference  in  computational  requirements  between 
N=15  and  N=16  was  the  first  indication  of  the  magnitude  of 
the  time  problem  which  would  be  encountered  in  dealing  with 
higher  values  of  N.  At.  N=17,  the  number  of 
cross-tcorrelation  pairs  increased  to  the  point  where  it  was 
also  necessary  to  re.vise  the  final  matching  algorithm  to 
prevent  memory  overflow,  compounding  the  time  problem.  With 
444  sequences  passing  the  autocorrelation  test,  a  gated  set 
of  length  272  was  obtained.  Several  sets  of  length  nine 
were  found  using  4+  hours  of  computer  time,  at  which  point 
program  execution  was  terminated.  It  is  estimated  that 
approximately  65  hours  would  be  required  for  a  complete 
search  with  the  algorithm  being  used. 

It  was  obvious  that  it  would  be  futile  to  continue  the 
search  tc  higher  values  of  N  using  the  same  procedure. 
Attempts  at  developing  a  more  efficient  algorithm  were 
unsuccessful.  The  32K  on-line  memory  of  the  SDS-9300  was  a 
serious  constraint,  but  even   without   it,   the   anticipated 
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execution  time  for  the  algorithms  envisioned  far  exceeded 
the  computer  time  that  could  be  made  available.  A  simple 
subterfuge  was  eventually  employed  to  obtain  partial  results 
for  higher  values  of  N. 

The  numter  of  sets  of  a  given  maximal  length  which  were 
being  found  for  higher  values  of  N  were  far  greater  than 
needed  in  order  to  offer  an  acceptable  degree  of  flexibility 
in  selecting  a  particular  set  for  use  in  an  operational 
system.  Many  of  the  sets  differed  only  in  the  replacement 
of  cne  or  two  unique  words,  and,  because  of  the  nature  of 
the  searching  process,  tended  to  be  uncovered  within  the 
same  time  frame.  tfhen  the  program  was  terminated  pricr  to 
completion  cf  a  search,  as  at  N=17,  not  only  was  it  unknown 
whether  or  not  sets  pf  greater  length  than  found  actually 
existed,  the  oijes  that  were  found  were  for  the  most  part 
different  combinations  of  the  same  slightly  larger  set  of 
unigue  words.  Given  the  fact  that  a  partial  search  must 
suffice,  it  seemed  appropriate  to  effect  a  reduction  in  the 
number  of  seguences  which  passed  the  autocorrelation  test  to 
the  point  that  the  cross-'correlation  search  could  be  run  to 
completion.  Doing  so  again  creates  uncertainty  as  to 
whether  or  not  sets  of  maximal  length  have  been  found. 
However,  given  the  excessive  number  of  sets  being  obtained, 
there  is  some  hope  that  some  of  them  might  still  appear  when 
the  autocorrelation  basis  is  reduced. 

Increasing  the  effective  detector  aperture  creates  a 
more  strict  autocorrelation  condition  and  thus  reduces  the 
numfcer  of  seguences  passing  the  autocorrelation  test.  A 
mathematically  artificial  increase  affects  detector 
performance  only  in  that  the  probability  of  false  detection 
is  reduced  (unnecessarily)  in  the  region  immediately  outside 
the  real  aperture.  The  effect  of  a  varying  aperture  was 
first  tested  using  N=15  and  N=16,  for  which  complete  results 
had  been  previously  obtained.  For  N=15,  the  maximal  length 
set  remained  at  three  for  aperture  widths  of  three  through 
six  bits.   The  number  of  such  sets  totaled  121  for   a   three 
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bit  aperture  versus  six  for  a  six  bit  aperture.  For  N=16, 
the  maximal  length  remained  at  six  for  aperture  widths  of 
three  through  six  bits,  with  corresponding  totals  of  724, 
345,  56  and  8.  Four  sets  of  length  five  were  obtained  for 
an  aperture  fcidth  of  seven  bits. 

Indeed,  up  to  a  point,  the  maximal  length  is  independent 
of  aperture  variation.  Partial  results  were  obtained  for 
values  cf  N  egualling  17,  18,  19  and  20  using  this 
technigue.  The  aperture  was  set  high  and  then  reduced  in 
increments  until  search  time  became  excessive.  The  results 
so  obtained  again  fell  short  of  expectations  in  that  much 
computer  time  was  st^.11  required  to  complete  the  partial 
search.  For  a  given  number  of  sequences  passing  the 
autocorrelation  test,  it  was  expected  that  search  time  would 
be  of  the  same  order  of  magnitude  for  the  different  values 
of  N.  In  fact,  search  time  increased  rapidly  with  N.  For 
example,  with  a  five  bit  aperture,  150  sequences  passed  the 
autocorrelation  test  at  N=16  and  the  total  search  time  was  3 
min,  50  sec.  With  a  16  bit  aperture,  142  sequences  passed 
the  autoccrrelation  test  at  N=20  and  the  search  was  only 
approximately  1/30th  complete  when  terminated  after  two 
hours.  The  imposition  of  a  stricter  autocorrelation 
criteria  produces  a  basis  of  the  same  order  as  for  a  lower 
value  of  N,  but  the  cross-correlation  properties  of  these 
seguences  is  such  that  many  more  combinations  of  match  sets 
exist.  That  is,  any  given  unique  word  cross-correlates  well 
with  many  more  seguences  in  the  N=20  basis  than  in  the  N=15 
basis. 

Data  on  the  numbers  of  sets  found  with  h=2  for  various 
values  of  N  and  various  aperture  widths  is  presented  in 
Table  4-2;  "W"  is  the  aperture  width  and  "NAUTO"  is  the 
numfcer  of  sequences  passing  the  autocorrelation  test.  An 
asterisk  denotes  partial  results  obtained  from  an  incomplete 
search.  Table  4-3  contains  similar  data  for  h= 1  and  h=3, 
and   serves   to  illustrate  the  effect  of  a  varying  threshold 
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on   the   numbers   involved.    A   representative  listing   of 
maximal  length  sets  .(for  h=2)  is  contained  in  Appendix  A. 


Gated  Set 

Set 

Number  of 

Maximal 

Maximal 

Maximal 

N 

a 

NAUTO 

Length 

Length 

iSM^ii  Se,ts. 

14 

4 

4 

0 

0 

- 

15 

i 

110 

2 

2 

2 

6 

26 

7 

3 

6 

5 

36 

13 

3 

16 

4 

60 

21 

3 

72 

3 

78 

22 

3 

121 

16 

1 

70 

20 

5 

4 

6 

110 

47 

6 

8 

5 

150 

66 

6 

56 

4 

190 

84 

6 

345 

3 

222 

95 

6 

724 

17 

9 

56 

21 

6 

2 

7 

260 

145 

9* 

unkn 

5 

444 

272 

9* 

unkn 

18 

13 

20 

14 

4 

18 

12 

46 

26 

6 

63 

11 

60 

47 

9 

4 

10 

142 

102 

10* 

unkn 

9 

358 

254 

unkn 

unkn 

19 

15 

36 

26 

5 

89 

14 

106 

83 

9* 

unkn 

13 

204 

157 

unkn 

unkn 

20 

16 

142 

113 

10* 

unkn 

15 

358 

306 

unkn 

unkn 

Table  4-2 
Sets  Obtained  For  h=2 
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10 

11 

12 
13 

19 

20 


w 

NA;0TO 

Gated  Set 
Maximal 
Length 

Set 

Maximal 
Length 

Numb€r  of 

Maximal 
Lengtjj  Sets 

7 
5 

4 
8 

2 
3 

2 

2 

5 
2 

9 
7 
5 

8 
14 
30 

2 

6 

13 

0 
3 
4 

10 
20 

9 
7 
5 

30 
44 
78 

14 

25 
44 

5 
6 

9 

16 

20 
10 

9 

7 
5 

32 
130 

182 

54 

84 

120 

10* 
11* 
unkn 

unkn 
unkn 
unkn 

5 

0 

0 

0 

- 

7 

5 

14 
7  6 

0 
17 

0 
4 

20 

Table  4-3 
Sets  Obtained  For  h=1  and  h=3 


D.   FALSE  DETECTION  PROBABILITIES 

The  probability  of  miss-detection  is,  fcr  constant  h, 
the  same  for  all  sequences  of  a  given  length. 
Representative  numbers  are  contained  in  Table  4-1.  Given  an 
acceptable  combination  of  length  and  detector  threshold,  it 
is  the  relative  values  of  false  detection  probability  that 
then  determine  whether  or  not  a  particular  sequence,  or  set 
of  sequences,  is  suitable  for  use  as  a  unique  word.  The 
auto-  and  cress-correlation  criteria  were  established  with 
this  in  mind.  As  a  point  of  final  analysis,  the  eguations 
of  Appendix  C  were  used  to  calculate  the  false  detection 
probabilities.  Representative  results  are  included  in  the 
following  paragraphs. 

1  •   Probability  of  Fjtlse  ^etget j,on  in  t h§   Random  R e ,g  ion 

In  an  efficient  TDMA  system,  the   number   of   random 

(information)   bits   in   a   frame   far  exceeds  the  number  of 

fixed  bits.   The  majority  of  terms  in  the  expression  fcr  the 

probability   of   one   or   more   false   detections  per  frame, 
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eguation  C.10,  are  therefore  of  the  form  given  by  equation 
C.4,  the  probability  of  false  detection  when  only  random 
bits  are  in  the  detector.   With  h=2, 

P{A(tH       =    3.693    x    10        for    N  =  15 

r 

-4 
=    2.012    x    10        for    N=20 

2-  Cross-Correlation   Baugd. 

When  cne  or  <more  unique  uord  bits,  other  than  the 
one  for  which  the  detector  is  set,  are  contained  in  the 
detector,  the  cross-correlation  criteria  establish  an  upper 
bound  on  the  probability  of  false  detection  (equations  3.4, 
3.6  and  3.7).  These  individual  event  probabilities  can  in 
turn  be  inserted  into  equation  3.S  to  obtain  an  upper  bound 
for  the  probability  pi  one  or  more  false  detections  (F1)  as 
any  one  unique  word  passes  through  the  detector.  Several 
such  calculations  were  performed;  for  N=15  and  N=20,  the 
results  were, 

-3 

F«     =  5.831  x  10    for  N=15 
max 

-3 
=  6.332  x  10    for  N=20 

Thus  F"  is  within  an  order  of  magnitude  of  the 
individual  event  probabilities  in  the  random  region. 
Because  of  the  t.ypically  large  number  of  random  events,  F* 
in  fact  makes  an  insignificant  contribution  to  the  total  F 
over  the  period  of  one  frame.  Given  the  fact  that  a 
cross-correlation  false  detection  pulse  is  potentially 
periodic,  an  "insignificant"  F'  is  precisely  what  is 
desired. 

3-  Exaffpj.e.    System   Calculations 

In  crder  to  compute  the  probability  of  one  or  more 
false  detections  over  a  frame,  it  is  necessary  to  have 
available   the   specific   frame  composition  of  the  operating 
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system  in  question.  A  system  was  hypothesized  for  the 
purpose  of  obtaining  a  general  feeling  for  the  relative 
magnitude  of  the  different  factors  involved.  The  number  of 
stations  was  set  at  three  in  order  to  simplify  the 
calculations  and  permit  selection  of  a  set  of  unique  words 
from  those  available  at  N=15.  The  following  set  was 
selected  at  random: 

(1)  000001100101111 

(2)  00010C100011111 

(3)  000101001101111 

This   set   possesses   the   following   correlation 
properties: 

AUTCCOREELATIQN  C^OSS-CSHSJ^ATION 

RJil 
t     JZlzill  J3J^JU  ill-Ill 

0  0 

0  0 

0  0 

0  1 

1  2 

2  2 

3  3 

4  2 

5  5 

4  5 

5  3 
8  6 
8  10 

6  8 
12  10 
10       10 

8  8 

8  8 

6  5 

3  6 

6  4 

3  5 

4      4 

3  3 

2  1 

1  1 

0  0 

0  0 

0  0 


The  hypothetical  system  is  patterned  after  the 
operational  system  described  in  Sef.  3.  The  transmission 
rate  is  60  Mbps  (4-piiase)  ,  the  frame  length  is  750  usee,  the 


£JU 

±1 

All 

S2X.  iJI 

0 

15 

15 

15 

1 

9 

9 

7 

2 

7 

7 

7 

3 

7 

7 

7 

4 

5 

6 

6 

5 

4 

4 

5 

6 

5 

3 

6 

7 

5 

5 

2 

8 

4 

4 

3 

9 

2 

2 

3 

10 

1 

1 

2 

11 

0 

1 

1 

12 

0 

0 

0 

13 

0 

0 

0 

14 

0 

0 

0 

-14 

0 

-13 

.  0 

-12 

0 

-11 

0 

-10 

0 

-9 

2 

-8 

4 

-7 

5 

-6 

5 

-5 

7 

-4 

6 

-3 

5 

-2 

7 

-1 

11 

0 

11 

+  1 

9 

+2 

9 

+  3 

9 

+  4 

3 

+  5 

3 

+6 

5 

+  7 

3 

+  8 

4 

+  9 

2 

+  10 

2 

+  11 

1 

+  12 

0 

+  13 

0 

+  14 

0 
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guard  band  is  24  bits  in  width  and  there  are  60  fixed  bits 
in  the  preamble  for  .carrier  and  bit  timing  recovery.  The 
latter  are  arbitrarily  assumed  to  consist  of  a  sequence  of 
30  1's  for  carrier  recovery,  followed  by  a  sequence  cf  30 
alternating  0's  and  1's  for  bit  timing  recovery.  Bits  are 
clocked  in  from  the  left,  with  the  detectors  set  tc  the 
unigue  wcrds  as  listed  above. 

It  is  first  necessary  to  determine  the  number  of 
random  versus  fixed  detection  events  per  frame.  The  product 
of  the  transmission  bit  rate  and  frame  length  gives  the 
number  of  bits  in  one  frame,  n=45000.  The  beginning  of  each 
station's  burst  is  cpmposed  of  75  fixed  value  bits  (recovery 
plus  unigue  word) .  When  one  or  more  of  these  bits  is  in  the 
detector,  the  detection  event  is,  by  definition,  "fixed." 
Therefore,  for  each  station's  burst,  there  are  75  fixed 
detection  events  as  each  bit,  from  the  first  to  the  last,  is 
decked  intc  the  detector;  plus  an  additional  15  fixed 
events  as  the  last  bit,  and  those  immediately  preceding  it, 
are  clocked  out  of  the  detector.  The  total  number  of  fixed 
detection  events  per  frame  is  thus  equal  to  270  (3x90)  .  The 
number  of  random  events  is  44730  (45000-270). 

With  this  number  of  random  events,  the  upper  bound 
on  the  probability  of  one  or  more  false  detections  per  frame 
in  the  random  region  alone  is  (equation  3.  10) , 

r3 
F  <  rP{A(t)}   =  44730(3.693  x  10   ) 

r 

F  <  165 

which  is  net  a  particularly  enlightening  result  for  a 
probability  upper  bound.  It  is  indicative,  however,  of  the 
fact  that  a  false  detection  is  almost  certain  during 
acquisition.   By  way  of  comparison,  the  fief.   11   expression 
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for  the  probability  of  no  false  detection  in  the  random 
region  gives, 

44730  -72 

1  -  F  =  [  1  -  P{A(t)}  J       =  1.3  x  10 

r 

Therefore,  it  can  be  concluded  that  many  false  detections 
will  occur  during  the  acquisition  phase,  and  that,  fcr  the 
system  design  hypothesized,  little  can  be  done  to  improve 
the  situation,  short  of  using  some  type  of  acguisition  gate. 
This  is  generally  true  for  the  typical  prototype  TDMA 
systems  which  have  been  developed.  No  attempt  was  made  to 
compute  the  expected  acguisition  time,  a  subject  treated 
extensively  in  Bef.  11. 

If  an  acquisition  gate  is  employed,  the  detector  is 
active  for  a  maximum  of  2B=48  detection  cycles  pricr  to  the 
receipt  of  the  associated  unigue  word.  All  of  these  are 
fixed  detection  events.  At  each  cycle,  the  autocorrelation 
values  can  be  added  to  the  number  of  matches  occurring  in 
the  recovery  bits  to  obtain  the  number  of  total  fixed 
matches,  k.  Equation  C.8e  can  then  be  used  to  obtain  an 
upper  bound  en  Fj   Dping  so,  for  detector  number  (2)  , 

-11 

F  <  3.996  x  10 

10 
1  >    2.503  x  10    frames  -  217  days 

For  detector  number  (3) , 

-4 
F  <  5.993  x  10 

T  >  1669  frames  =  1.25  sec 

The  mean  time  between  false  detections  (T)  is  more 
than  acceptable  for  both,  indicating  a  high  probability  of 
achieving  acquisition  in  just  one  pass  of  the  unigue  word. 
Acquisition  cf  the  reference   unique   word,   assumed   to   be 
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number  (1)  here,  must  still  be  effected  without  benefit  of 
gating,  and  thus  remains  a  potentially  time  consuming 
process. 


V.   CONCLUSIONS 

This  research  effort  was  undertaken  for  the  express 
purpose  of  finding  unique  word  seguences  which  would  be 
suitable  for  use  in  a  TDMA  satellite  communication  system. 
The  correlation  properties  of  potential  unigue  words  were  of 
primary  concern  in  this  effort,  criteria  being  established 
which  would  insure  a  low  probability  of  false  detection. 
Numerous  sets  of  unique  words  were  found  using  these 
criteria. 

The  results  obtained  can  not,  unfortunately,  be  applied 
to  any  given  operatipnal  TDMA  system  without  some  additional 
consideration  of  false  detection  probabilities.  In  an 
attempt  to  generalize  findings,  it  was  necessary  to  ignore 
the  fact  that  each  system  has  its  own  peculiar  set  of  fixed 
value  bits,  in  addition  to  the  unigue  words,  which  are 
transmitted  in  each  station's  burst.  These  fixed  value  bits 
are  in  fact  an  extension  of  the  unique  word  sequence  and 
produce  a  cross-correlation  distribution  of  their  oi»n  when 
passing  through  a  detector.  The  cross-correlation  may  very 
well  be  such  that  a  set  of  unigue  words  selected  at  random 
from  this  work  would  prove  unacceptable.  Fortunately,  the 
number  of  sets  available  of  any  given  length,  together  with 
the  existence  of  the  complementary  options,  provides  a 
flexibility  sufficie&t  to  cover  any  foreseeable  situation. 

The  length  of  the  sets  found  fell  short  of  expectations. 
To  begin  with,  in  setting  cross-correlation  criteria  based 
on  worst-case  possibilities  for  the  number  of  detector 
agreements,  unique  words  which  might  actually  be  useable  in 
a  particular  system  were  excluded.  That  is,  the  number  of 
agreements  produced  when  non-unique  word  fixed  bits  are 
stepping  through  the  detector  can  seldom  be  expected  to 
approach  the  worst-case  limit.  This  fact  was  noted  at  the 
outset,   however,   and   accepted   in  order  to  generali2e  the 


47 


results  as  much  as  possible.  The  real  disapp ointing  factor 
was  that,  even  with  the  strict  correlation  criteria,  a 
computer  search  could  not  be  designed  to  uncover  with 
certainty  the  sets  of  maximal  length  for  values  of  N  in 
excess  of  sixteen. 

There  is  a  need  for  an  efficient  algorithm  for  comtining 
unigue  wcrds  into  sets  possessing  good  cross-correlation 
properties.  As  stated  in  Ref.  11,  the  problem  is  a  pu2zling 
one..  It  is  a  problea  which  can  undoubtedly  be  solved  only 
by  adopting  some  type  of  short-cut  technigue,  wherein  but  a 
few  of  all  the  multitude  of  different  matching  combinations 
are  considered.  The  ability  to  predict  the  maximal  length 
for  a  given  N  would  be  of  tremendous  benefit  in  this  regard. 

The  ultimate  achievement  would  be  the  discovery  of  an 
underlying  structure  which  enabled  one  to  extract  unigue 
word  seguences  without  the  need  to  resort  to  brute  force 
computer  search  technigues. 
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APPENDIX  A 
UHI2UI!  W2RJ2  JETS 

Selected  sets  pf  unique  words,  representative  cf  the 
many  found,  are  presented  herein.  The  listing  for  N=15  is 
complete.  In  the  interest  of  space,  unique  words  are  listed 
in  octal  fcrmat.  All  sets  are  valid  for  a  detector 
threshold  setting  of  two. 

A  "gated  set"  is  one  in  which  the  first  sequence  listed 
(the  supposed  reference)  cross-correlates  well  with  all  the 
sequences  which  follow.  A  "complete  set"  is  one  in  which 
each  sequence  cross-correlates  well  with  every  other 
sequence. 


N  =  15 


MAXIMAL  1ENGIH  GATED  SET 


01457   01347   02067   02327   02367   02467   03067 
04437   05157   04157   04657   04517   04117 


COMPLETE  SETS 


(1) 

00627 

02067 

02217 

(?) 

06627 

022^7 

02327 

(3) 

06627 

02217 

02367 

(<0 

01317 

01667 

02327 

(5) 

01457 

04437 

05157 

(6) 

01457 

04437 

04157 

(7) 

01457 

04437 

04117 

(8) 

01667 

02327 

02637 

(9) 

01667 

02367 

02717 

(10) 

01667 

02367 

02637 

(11) 

01667 

03367 

02637 

(12) 

05477 

04757 

03557 

(13) 

05477 

03557 

05157 

(14) 

05477 

03557 

C4157 

(15) 

03137 

04437 

05157 

(16) 

04437 

04157 

C3057 
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MAXIMAL  LENGTH  GATED  SET 


N=16 


002637 


001627 
002367 
003327 
004167 
004657 
006467 
014337 
011637 
010757 
010157 


001647 
002467 
003347 
004327 
004717 
006567 
010337 
011237 
012357 
011657 


01667 
02647 
003367 
004367 
004727 
006767 
011537 
014237 
010357 
006657 


002247 
003067 
003557 
004517 
005167 
014277 
014537 
010237 
006557 
010657 


002327 
003127 
003667 
004617 
006167 
011737 
006137 
006757 
005157 
006257 


010457   010717   011517   010517   010617 


002347 
003267 
004157 
004627 
006267 
012337 
014137 
004757 
006157 
01  1457 


COMPLETE  SETS 


(1) 

001457 

002717 

005317 

011237 

010657 

011517 

(2) 

001667 

002637 

004657 

005167 

006467 

006157 

(3) 

002367 

005167 

006467 

005477 

0  06137 

006257 

(<0 

002467 

002637 

003267 

004367 

006137 

006257 

(5) 

002637 

003267 

004367 

006467 

006137 

006257 

(6) 

002637 

003267 

005167 

006467 

006137 

006257 

(7) 

C02637 

004657 

006137 

01123  7 

010357 

011457 

(8) 

002717 

005317 

010337 

003137 

010657 

011517 

(9) 

003057 

004567 

005067 

006267 

004437 

005157 

(10) 

003267 

004367 

004717 

005467 

003137 

005157 

(11) 

003557 

011537 

005637 

010557 

011657 

011317 

(12) 

004567 

OQ5067 

006267 

006137 

004437 

005157 

(13) 

004417 

003.137 

004437 

005157 

010157 

011517 

(14) 

005317 

011077 

010337 

003137 

010657 

011517 
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N  =1  7 


COMPLETE    SETS 


(1) 

001667 
014457 

006257 
012157 

011057 
011437 

011167 
014237 

012467 

(2) 

003137 
011767 

005567 
012667 

005717 
012717 

007157 
013167 

007327 

(3) 

003327 
012667 

005567 
012717 

007157 
013167 

007267 
021017 

011767 

(4) 

002637 
024.157 

021137 

022257 

014137 
014257 

011437 
021457 

024437 

(5) 

0027  1  7 
007317 

003237 
011727 

005567 
012267 

006327 
013167 

007267 

(6) 

.003137 
011767 

005667 
012267 

005717 
012717 

006727 
013167 

007157 

C?) 

005567 
01172  7 

005717 
012667 

006657 
012717 

007157 
013167 

007267 

(8) 

005717 
011767 

006657 
012667 

006727 
012717 

007267 
013167 

007567 

N  =  1fi 


COMPLETE  SETS 


(1) 


(2) 


(3) 


M 


(5) 


00145  7 
020567 

00465  7 
016667 

012717 

006277 

014657 
053637 

021557 
050637 


004767 
021667 

004767 
021767 

014657 
014337 

030277 
045437 

044177 

045437 


01  1637 
025517 

005717 
026317 

025517 

045437 

042477 
042757 

042677 
041357 


014657 
014337 

012717 
006277 

044177 
041357 

043337 
022357 

024677 
022357 


015527 
024637 

015727 
014337 

040677 
022357 

050337 
044357 

030537 
031257 
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COMPLETE  SETS 


N  =  19 


(1) 
(2) 


(3) 


CO 


0004657 
0010767 

0004657 
0104577 

0007557 
0104577 

0Q07557 
0104577 


0005567 
0011167 

0005637 
0.106677 

0011537 
0106677 

0013357 
0101177 


0005637 
0030577 

0007557 
0014277 

0013357 
0014277 

0013477 
0106677 


0006767 
0042277 

0013357 
0101477 

0013477 
0101477 

0015757 
0014277 


0007557 
0022477 

0015757 
0022477 

0015757 

0022477 

0046577 
0101477 


COMPLETE  SETS 


N=2Q 


(1) 
(2) 


(3) 


(4) 


(5) 


0004437 
0061177 

000<*437 
0021137 

0007137 
0032737 

0013237 
0035737 

0021137 
0064377 


0006277 
0054677 

0007137 
0024337 

0011077 
0060577 

0013477 

0062577 

0023077 
0035577 


0023477 
0024677 

0011537 
0C3  0737 

0013237 
0054677 

0014277 
0050677 

0024477 
0054577 


0030737 

0C43277 


0030277 

0015637 
0047277 

0020677 
0042277 

0024737 
0036677 


0064577 
0040477 

637 
0045477 

0021337 
0041477 

0032737 
0047477 

0030337 
0043477 
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£PPEi22Ix  b 

flffTEgMINATION  OF  BJRST  DRIFT  JgATE 

In  the  absence  of  transmitter  burst  position  corrections 
and  under  zero  doppler  conditions,  a  given  unique  word  will 
be  received  every  T  seconds.  (Other  factors,  such  as  clock 
stability  and  logic  jitter,  affect  this  periodicity; 
however,  they  are  inconsequential  compared  to  potential 
doppler  effects  and  will  be  ignored.)  Burst  drift  is 
defined  as  the  doppler  induced,  real-time  shift  of  a  burst 
from  its  periodic  position. 

Burst  drift  rate  is  a  function  of  the  relative  velocity 
between  transmitting  station,  satellite  and  receiving 
station.  Interest  here  is  in  determination  cf  the  maximum 
drift  rate  that  can  be  expected.  First,  let  us  consider  the 
motion   of   the   typical   communication   satellite    in   an 

remains  at  the  same  altitude  over  a  given  point  on  the 
equator.  In  reality.,  it  is  physically  impossible  to  achieve 
and  maintain  a  perfectly  synchronous  orbit  [12].  Reference 
13  derives  equations  which  can  be  used  to  estimate  the 
worst-case  relative  velocity  between  a  geosynchionous 
satellite  and  a  given  point  on  the  earth's  surface.  The 
orb^t  eccentricity  and  inclination  are  the  two  primary 
factors  which  influence  the  satellite's  motion. 

A  non-circular  orbit  results  in  a  sinusoidal  variation 
of  the  satellite  altitude.  To  a  first  order  approximation, 
the  resulting  doppler  shift  is  independent  of  earth  station 
coordinates.  If  e  is  the  eccentricity  of  the  orbit,  the 
maximum  relative  velocity  is  approximately, 


v   =  12100e  Km/hr  =  3361e  m/sec 
e 


53 


Using  .0004  frpm  Table  1  of  Ref.  14  as  a  typical 
eccentricity  for  later  generation  communication  satellites, 

v   =1.3  m/sec 
e 

When  the  orbit  is  inclined  from  the  eguatorial  plane, 
the  satellite  moves  in  a  figure-8  pattern,  elongated  north 
and  south.  The  north-south  motion  predominates  and  is 
sinusoidal  in  nature.  If  i  is  the  orbit  inclination  in 
radians,  the  resulting  maximum  relative  velocity  as  seen  by 
a  station  at  latitude  L  can  be  approximated  by, 

v      =   2092^sin(L)     Kia/hr    =    581isin(L)    m/sec 
i 

Satellite   visibility    extends   to   approximately   80° 

latitude.   At  this   latitude,   the   satellite   appears   just 

above   the   horizon   and  relative   velocity   due   to   orbit 

inclination  is  a  maximum.   (It  is  in   fact   highly   unlikely 

that  a  station  would  attempt  to  use  a  satellite  below  10°  in 

elevation.)   Thus,  assuming   a   maximum   inclination   cf   3° 

(.052  rad)  and  L=80°, 

v   =  29.7  m/sec 
i 

A  pessimistic  estimate  of  the  maximum  relative  velocity 
due  to  satellite  motion  alone  can  be  obtained  by  combining 
the  effects  of  orbit  eccentricity  and  inclination,  giving, 

.v   =  v   +v   =31  m/sec 
s    e    i 

Next,  the  contribution  of  earth  station  motion  must  be 
considered.  Reference  13  analyzes  the  headway,  roll  and 
pitch  components  of  ship  motion  and  concludes  that  the 
maximum  possible  relative  velocity  contribution  is, 

v   =38  m/sec 

g 
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For  a  shipboard  platform, then,  the  maximum  relative 
velocity  between  ground  station  and  satellite  is, 

v=v   ♦  v   =69  m/sec 
g    s 

Kith  representative  numbers  for  satellite  and  surface 
platform  motion  in  hand,  it  is  next  necessary  to  consider 
how  they  night  combine  to  produce  a  shift  in  the  position  of 
a  received  unigue  word.  The  case  of  a  station  receiving  its 
own  unigue  word  (transmitted  up  to  the  satellite  and  back) 
is  a  simple  two-station  relative  motion  problem  and  the 
well-kncwn  dcppler  eguation  applies. 

d  =  2vf/c  (B.5) 

where   f   is   the   transmitted   bit   rate  and  c  the  speed  of 
light. 

The  case  of  a  station  receiving  another's  unigue  wcrd  is 
a  three-station  motion  problem.  Maximum  doppler  occurs  when 
the  maximum  relative  velocity  exists  between  both  ground 
stations  and  the  satellite  simultaneously;  that  is,  when 
maximum  up-link  doppler  is  combined  with  maximum  down-link 
doppler.  Such  a  situation  is  produced  when  both  ground 
stations  are  at  the  same  80°  latitude  and  both  are  either 
opening  or  closing  the  satellite  at  the  maximum  rate.  For 
this  case  also,  then,  with  maximum  relative  velocity,  v, 
existing  between  both  stations  and  the  satellite,  eguation 
B.5  applies. 
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APPENIilX  C 
DETECTION  PROBABILITIES 

Expressions  for  the  probability  of  missed  and  false 
detections  are  derived  in  the  following  paragraphs.  These 
expressions  depend  somewhat  on  general  detector  operation, 
but  not  to  the  extent  that  their  application  is  restricted 
narrowly  to  any  one  particular  system. 

A.   MISS-EETICTION  PROBABILITY 

If  the  bit  error  rate  (BEE),  i.e.  the  probability  that  a 

received   bit   will  be  in  error,  is  p  and  the  number  of  bits 

in   a   unigue   word   is   N#   then,   assuming   errors   occur 

independently,   the   probability   of   a   unigue   word   being 

N 
received  error-free  4-s  (1-p)  .   With  the  detector   threshold 

set  at  h,  a  unigue  word  will  be  detected  if  it  contains  h  or 
fewer  errors.  The  probability  of  detection,  P,  is  thus 
obtained  by  summing  terms  in  the  binomial  distribution  which 
correspond  to  the  prpbability  of  no  errors,  the  probatility 
of  one  error,  up  to  the  probability  of  h  errors.  Using 
C(N,i)  to  represent  the  binomial  coefficient  N!/i;(N-i);,  P 
is  given  by, 

h        i     N-i 
P  =.2n    C(N,i)p  (1-p)  (C.1) 

1=0 

The  miss-detection  probability  (the  probability  of 
receiving  a  unigue  word  with  more  than  h   errors)   is   given 

by,  . 

N  i      N-i 

Q    =  1-P  =  2  C(N,i)p  (1-p)         (C.2) 

i  =  h+1 

The  mean  time  between  miss-detections  is  egual  to  the 
reciprocal  of  the  expected  value  of  "number  of  misses  per 
unit   tice."    As   there  is  only  one  possible  miss-detection 
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event  per  frame,  .the  expected  value  is  equal  tc  the 
miss-detecticn  probability,  Q,  and  the  mean  time  between 
miss-detecticns  is, 

"1  ~1 

a   =  Q   frames  =  TQ    sees        (C.3) 

B.   PfiOBAEILITY  OF  FALSE  DETECTION 

The   following   notation   will   be   used   in   developing 

expressions  for  false  detection  probabilities: 

A  (t)       A  false  detection  event  at  time  t. 

Used  interchangeably  with  "A  .  •• 

t 

F{A(t)}    The  probability  of  A  (t)  . 

E  (t)       No  false  detection  at  time  t. 

Used  interchangeably  with  "B  ." 

t 

F{B(t)}    The  probability  of  B  (t)  . 

t         An  integral  time  variable,  the  basic  unit 

of  which  is  one  bit  cycle  (t=ib) . 


On  any  given  bit  cycle,  other  than  the  one  in  which  the 
actual  unigue  word  is  completely  contained  in  the  detector, 
a  false  detection  pulse  will  be  produced  if  N-h  or  mere  of 
the  N  bits  in  an  active  detector  match  with  the  detector 
setting  (B>h) .  The  probability  of  this  occurring  depends  on 
whether  cr  net  the  value  of  transmitted  bits  are  fixed  or 
variable.  Ihe  beginning  portion  of  each  burst  is  composed 
of  fixed-value  bits  (timing  recovery  and  unigue  word) ,  while 
the  remainder  of  the  burst  contains  bits  of  variable  value. 
The  probability  of  one  of  the  latter  being  transmitted  as  a 
1  01  or  a  '1'  depends  on  the  type  of  information  transmitted. 
As  a  first  order  approximation,  it  will  be  assumed  that  the 
variable  tits  are  in  fact  random,  taking  on  the  value  ,0I 
and  M1  with  egual  pxobability .   With  random  bits   alone   in 
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the      detector,      the      probability  of  a    false   detection    en   any 
given    bit   cycle   is, 

P{A(t)}      =.J      C(N,i)  (1/2)l(1/2) 

r      i=0 

N      h 
=     (V2)       ?      C(N,i)  (C.4) 

1=0 

The  probability  of  false  detection  is  not  so  easily 
computed  when  there  are  fixed  value  bits  in  the  detector. 
"Fixed  value"  is  used  to  denote  those  bits  in  the  preamble 
which  are  transmitted  with  the  same  value  {'0'  or  *1')  in 
each  frame.  Because  of  the  presence  of  transmission  errors, 
these  bits  may  actually  be  received  with  the  opposite  value, 
the  probability  of  which  is  p.  An  expression  fcr  the 
probability  of  false  detection  can  be  derived  based  on  the 
number  of  f4.xed  bits  in  the  detector  and  the  number  that 
matched  with  the  detector  setting  when  transmitted.  Such  an 
expressicn  must  account  for  all  the  possible  combinations  cf 
errors  within  the  fixed  bits  and  as  such  will  include 
additive  terms  containing  higher  powers  of  p.  With  p<<1,  a 
good  approximation  can  be  obtained  by  discarding  all  terms 
which  contain  other  than  the  lowest  power  of  p.  Therefore, 
rather  than  derive  a  complex  general  expression  from  which 
an  approximation  can  be  drawn,  the  condition  p<<1  will  be 
invoked  at  the  outset  in  order  to  simplify  the  derivation. 

Suppose  there  are  m  fixed  bits  in  the  detector  and  that 
k  cf  them  matched  with  the  detector  setting  when 
transmitted.  A  false  detection  will  occur  when  the  number 
of  agreements  within  the  fixed  bits,  plus  the  number  of 
agreements  within  the  remaining  N-m  random  bits  is  greater 
than  or  egual  to  N-h,.  With  p<<1,  interest  can  be  focused  on 
those  particular  cpmbinations  which  reguire  the  fewest 
number  cf  transmission  errors  in  order  to  produce  a  false 
detection.   Three  separate  cases  will  be  considered. 
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CASE  1.  The  number  of  fixed  agreements,  kr  is 
greater  than  or  equal  to  N-h.  A  false  detection  is  assured 
if  the  m  fixed  bits  are  received  error  free,  and  the 
probability  cf  false  detection  can  be  approximated  by, 

P-{Mt)}   =  (1-P)  ;    k  +  h>N      (C.6a) 

CASE  2.  The  number  of  fixed  agreements,  k,  is  less 
than  N-h,  but  enough  random  bits  (N-m  of  them)  remain  to 
produce  a  total  of  N*h  or  wore  agreements.   That  is, 

k+N-m>N-h 

or, 

k  *  h  >  m 

A  false  detection  will  be  produced  when  there  are  no  errors 
in  the  fixed  bits  and  N-k-h  or  more  agreements  withic  the 
random  bits;  or,  eguivalently ,  when  there  are 

N-m-  (N-k-h)  =k  +  h-m 

or  fewer  disagreements  within  the  random  bits.  The 
probability  of  false  detection  can  thus  be  approximated  by 
multiplying  the  probability  of  receiving  the  m  fixed  bits 
error  free  by  the  probability  that  the  N-m  random  bits  will 
contain  k+h-m  or  fewer  disagreements. 


P{A(t)}   = 

f 

m      N-m   k+h-m 
(1-p)  (1/2)     2  C(N-m,i);    m<k+h<N     (C.Eb) 

i=0 


CASE  3.  The  number  of  fixed  agreements,  k,  plus  the 
number  cf  random  bits,  N-m,  is  less  than  N-h;  i.e.,  k«-h<m. 
A  false  detection  can  not  occur  unless  there  are  errors  in 
the  fixed  bits.  The  fewest  number  of  errors  are  reguired 
when  all  cf  the  random  bits  and  all  of  the  original   k   bits 
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match  with  the  detector  setting.  In  this  case,  enough 
errors  must  occur  in  the  remaining  m-k  fixed  bits  to  produce 
N-h  agreements  (the  minimum  required  for  false  detection) . 
The  number  of  errors  required  is, 

a  =  N  -  h  -  (k  +  N-m) 

=  m  -  k  -  h 

The  protability  of  false  detection  can  thus  be  approximated 
by  taking  the  product  of  (1)  the  probability  that  the  N-m 
random  tits  will  all  match  with  the  detector  setting,  (2) 
the  probability  of  receiving  the  original  k  matching  bits 
error  free  and  (3)  the  probability  that  the  m-k  fixed  tits 
will  contain  ,a'  errors. 


a-m     k         a     ro-k-a 
P{A(t)}   =  (1/2)    (1-p)  C(m-k,a)p  (1-p) 

N-m         a     ra-a 
=  (1/2)    C(m-k,a)p  (1-p) 

N-m  m-k-h     k+h 

=  (1/2)    C  (m-k,  m-k-h)  p      (1-p)    I         k«-h<m 

(C.8c) 


If  ail  tiie  bits  in  the  detector  are  fixed,  the 
probability  of  false  detection  in  any  given  bit  cycle  is 
given  by  eguatior^s  C.8a  and  C.8c  with  in=N. 


N 
P{A(t)}   =  (1-p)  ;    k+h>N       (C.8d) 

N-k-h     k+h 
P{A(t)}   =  C(N-k,  N-k-h)  p      (1-p)    ;    k+h<N  (C8e) 


1 •   ^oi » t  Probabilities 

Zguations  C. 4  and  C.8  give  the  probability  of  false 
detection  for  a  single  detection  event  within  the  frame. 
Calculation   of   the   probability   of   false   detection  over 
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several  consecutive  bit  periods,  specifically  over  one 
frame,  is  complicated  by  the  fact  that  successive  detection 
events  are  net  independent;  i.e., 

P{A(t),A(t+1)}*P{A(t)}P{A(t  +  1)}. 

That   this   is  true   can   be  seen  by  noting  that  successive 

detection  events  are  generated  by  shifting  the  contents   of 

the   shift   register   by  one  bit.   If  a  false  detection  does 

(or  does  not)  occur   at   one   step,   possibilities   for   the 

contents  or  the  shift,  register  are  necessarily  restricted  to 

certain  combinations  of  bit  values.   These  bit  values   carry 

over  to  the  next  detection  event  and  reduce  the  sample  space 

accordingly.   As  a  sample,   extreme   example,   consider   the 

case   of   the   unique  word  101010.   The  probability  of  false 

detection  with  random   bits   in   the   detector   and   a   zero 

6 
threshold   setting   is,  from  equation  C.U,  (0.5)  .   However, 

given  that  a  false  detection  occurs  (the   detector   contains 

101010),   the   detector  will  contain  X10101  en  the  next  step 

and  the  probability  of  false  detection  is  zero. 

Dependent  joint  probabilities  are  best  dealt  with  by 
locking  for,  and  finding,  grounds  for  the  assumption  of 
independence.  Schrempp  and  Sekimoto  [11]  based  their 
independence  assumption  on  the  previous  work  of  Hill  [1^]- 
The  validity  of  this  assumption  is  nearly  impossible  to 
access  without  making  lengthy  comparative  calculations  for 
unigue  words  of  different  length,  a  task  of  monunental 
proportions.  Rather  than  attempt  to  justify  such  an 
assumption  here,  jpint  probabilities  will  be  expressed  in 
terms  of  an  upper  bound. 

Ihere  are  n  detection  events  per  frame  during 
acquisition,  n-1  of  which  may  produce  false  detections. 
Defining  t=0  at  the  point  when  the  associated  unique  word  is 
completely  contained  in  the  detector,  the  probability  cf   no 
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false  detection  during  a  frame  period  is  given  by, 


H  =  P{E  ,B  ,-  —  ,B   }  (3.9) 

1   2      n-1 

=    P{£    }P{3    /B    }P[B    /B    B    }«««P{B         /B    B    •••B         } 
1  2      1  3       12  n-1       1    2  n-2 


Cr,      from      a      different      viewpoint,      if      F      is        the 
probability    cf   one    or    more    false   detections   per  frame, 

H    =    1    -   F 


n-1  n-1      n-1 

*   -.2,         P(A.)    -.2,       .2.         P(A.,A.) 
1=1  l         1=1       3=1  l      3 

n-1      n-1      n-1 
♦.2.        2.       .2.        P(A.,A.,A    )    - 
1=1       3=1      k=j  13k 

±    P(A    ,A    , »«»,A         ) 

12  n-1 


For  straight-forward  calculations,  F  is  obviously 
orders  of  magnitude  more  difficult  to  obtain  than  is  ft. 
However,  its  form  does  offer  possibilities  for  reducing  the 
complexity  cf  .the  problem.  Although  successive  events  are 
not  independent,  events  separated  by  more  than  N  detection 
cycles  are  -  -  after  N  detection  cycles,  all  bits  frcm  the 
first  event  have  been  shifted  out  of  the  detector.  For 
j>i+N, 

P{A  ,A  }  =  P{A  }P{A  } 
i   J        i     3 

and  similarly  for  higher  order  joint  probabilities.  If 
P{A(i))<<1,  these  product  terms  are  small  compared  to  the 
individual  terms  from  the  first  summation  in  F.  However, 
for  the  dependent  .joint  probabilities  it  is  quite  possible 
that 

P{A  ,A  }  =  P{A  }P{A  /A  } 
i   j        1     j   1 
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is   of   the   same   order  as  P  {A  } .   The  same  can  be  said  for 

i 

higher  crder  terms,  but  with  a  decreasing  expectation  cf  its 

actually  prcving  true.  It  is  therefore  not  possible  to 
discard  as  insignificant  all  of  the  terms  involving  joint 
probabilities.   The  most  that  can  be  said  is  that, 

P{A  1>P(A  ,A  }>P{A  ,A  ,A  }  •  •  • 
1     i  3  i  3      k 

And  that, 

n-1 
1=1       l 

The  false  detection  probabilities  P{A(i)}  are 
obtainable  from  equations  C.U  and  C.8  as  appropriate.  For 
the  random  detection  events  (eguation  C.4) , 
P {A  (i)}=t  {A  ( j)} ;  for  the  fixed  detection  events,  the 
probabilities  must  be  computed  individually.  If  there  are  r 
random  false  detection  events  and  f=r.-'!-r  fixed  false 
detection  events  per  frame, 

F  =?     P{A.} 

max   i=1      i 

=  rP{A(t)j   «-  ^P{A(t)}        (C.10) 
r    r       f 

When  the  detector  is  gated  during  acquisition,  as 
described  in  Section  II. C,  the  detector  is  active  for  only 
2B+1  bit  cycles  about  the  nominal  position  of  the  associated 
unigue  wcrd.  The  true  detection  pulse  may  occur  anywhere 
within  the  aperture,  depending  upon  the  amount  of 
displacement  of  the  unique  word  from  its  nominal  position, 
and  the  maximum  number  of  false  detection  events  is  equal  to 
2B.  The  probability  of  one  or  more  false  detections  in  the 
aperture.  prior  tc  the  receipt  of  the  unique  word  is 
therefore  of  the  same  form  as  equation  C.10,  with  the  number 
of  events  egual  to  2B  vice  n-1. 
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In  the  steady-state  mode,  the  aperture  width  is  2D+1 
and  the  maximum  number  of  false  detection  events  is  equal  to 
2D. 

In  summary,  equation  C.10  provides  a  bound  on  the 
probability  of  one  or  more  false  detections  over  one  frame 
period,  or,  with  a  chanqe  of  limits,  over  an  acquisition  or 
steady-state  detection  aperture.  This  bound  is  useful  in 
estimating  the  frequency  of  false  detections  and  in 
assessing  the  relative  merits  of  different  unique  word 
comtinaticns . 

As  an  aside,  it  is  interesting  to  compare  the  random 
event  term  of  equation  C.10  with  the  results  of  Ref.  VI, 
obtained  by  assuming  independent  events.  For  randcm  bit 
detection  events,  it  was  concluded  that, 


F  =  1  -  [1  -  P{A(t)}  ]r 


Expanding, 


2 
I-  =  1  -  [1  -  rP{A(t)}  +  C(r,2)  P  {A  (t) }   -  •••] 


For  I{A(t)}«1,  a  basic  assumption  cf  Ref.  14,  an 
approximate  expression  for  F  is, 

F  =  rP{A  (t)} 

Which  is  identical  to  the  random  bit  term  in  equation   C.10. 
2.   Mean  Time  ^etyeei).  False  Detections 

A  false  detection,  both  in  the  acquisition  and 
steady-estate  modes,  deactivates  the  detector  and  sets  an 
aperture  for  the  next  frame.  The  probability  of  cycling 
through  several  successive  detection  events  without  arry. 
false  detections  is  therefore  of  primary  ccncern.  Fron  the 
opposite  viewpoint,  it  is  the  probability  of  one  or  more 
false  detections  over  a  given  time  span  which  is  meaningful, 
not  the  probability  of  false  detection  at  any  given  step. 
This  distinction  is  important  in  deciding  upon  what   tc   use 
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as   an   appropriate   measure   of   "mean   time   between  false 
detections"  (MTBFD) ,  and  in  turn  interpreting  its  meaning. 

If  the  detector  were  not  deactivated  by  a  false 
detection  pulse,  several  false  detections  might  occur  over 
the  pericd  of  one  frame.  The  MTBFD  would  then  be  equal  to 
the  reciprocal  of  tke  expected  value  for  "number  of  false 
detections  per  frame/."  With  detector  deactivation,  only  one 
false  detection  per  frame  can  actually  occur  and  the  MTBFD 
is  affected  accordingly.  In  this  case,  for  computational 
purposes,  the  probability  of  a  real,  physical  false 
detection  pulse  occurring  somewhere  within  the  frame  is 
egual  to  F,  the  probability  of  one  or  more  false  detections 
per  frame.  The  MTBFD  is  therefore  egual  to  the  reciprocal 
of  F  and  can  be  bounded  by  substitution  from  equation   C.10. 


-1     ^1 

3?  =  F    >  F    frames  (C.  11) 

max 
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A£PENDI£  J2 
COMPUTE  f ROGRAjJ 

The  computer  program  was  written,  and  debugged,  in  three 
separate  functional  segments  and  then  patched  together  to 
form  one  continuous  program.  These  segments  will  be 
described  and  listed  separately.  The  last  segment  listed  is 
the  basic  cross-correlation  matching  routine.  Several 
different  revisions  were  employed  in  an  attempt  to  overcome 
time  and  memcry  limitations. 


Arrays 


GUW(IEIM)  -  Seguences  which  pass  the   autocorrelation 
test. 

GCTUW(IDIM)  -  Complement  transpose  of  GUW. 

GTUW(IDIM)  -  Transpose  of  GUW. 

GCUW  (IDIM)  -  Complement  of  GUW. 

MTEME(IDIM2)  -  Used  for  temporary  storage  of  good 
cress-correlation  matches. 

INCE  (    )  -  Used  to  store  pointers  needed  in  the 

cress-correlation  matching  process.   Dimension  greater 
than  or  equal  io  maximal  length  set  expected. 

MATCH  (IDIK1,IDIM2)  -  Matrix  in  which  good 
crcss-*correlat:j.on  matches  are  recorded. 

MCEOSS  (2IDIM)  -  Seguences  from  autocorrelation  sort 
which  are  checked  for  cross-correlation  (the  GUW  and 
GCTDH  subsets),. 

Equivalence  Statement  -  The  first  half  of  MCROSS  is 
stored  in  same  memory  location  as  GUW(  );  the  second 
half  in  same  location  as  GCTUW  (  ). 
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Variables  Common  to  All  Segments 

UW  -  Seguence  under  test. 

SUH  -  Shifted  version  of  UW.   A  seguence  is  shifted  one 
bit  at  a  time  and,  at  each  step,  added  mod-2  to  UW. 

TUW  -  Transpose  of  UW. 

CTUW  -  Complement  transpose  of  UW. 

CUW  -  Conplement  of  UW. 

N  -  Length  of  sequences  being  tested  (Input) . 

IEPS  -  Thresholds  h  of  thesis  main  body  (Input) . 

IALFA  -  Aperture  parameter;  equals  2D  (Input). 

NUW  -  Number  of  good  UW*s  obtained  in  autocorrelation 
sort. 

NAUTO  -  Number  of  UN's  and  CTUW's,  combined,  which  are 
checked  for  crpss-correlation;  equals  2NUW. 
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-  Autocorrelation  Sojrt 
Variables 

K  -  Smallest  integer  equal  to  or  greater  than  N/2. 

LAST,IAS11  -  Intermediate  skip  points  in  the  sequential 
checking  process.   Sequences  which  must  be  checked 
tend  to  occur  in  blocked  subsets,  with  two  redundant 
seguences  separating  each. 

IZ1  -  Next  to  last  sequence  which  must  be  checked. 

IZ2  -  Last  sequence  which  oust  be  checked  to  insure  all 
possible  non-redundant  sequences  have  been  considered. 

MASK  -  Binary  constant  consisting  of  N  1's.  right 

justified.   Used  to  obtain  the  complement,  of  an  N  bit 
seguence. 

IF1AG  -  Set  to  one  when  last  sequence  being  checked. 
Keys  the  termination  of  autocorrelation  search. 

JFLAG  *  Set  to  one  when  autocorrelation  condition 
reguires  R(t)=Q.   Autocorrelation  computation  for  a 
particular  UH  is  complete  when  JFLAG=IEPS. 

I  -  Shift  variable.   Equivalent  to  time  variable  of 
thesis  main  body. 

ICOUNT  -  The  autocorrelation  value  at  step  I. 
Output 

The  sequences  which  pass  the  autocorrelation  test,  and 
the  corresponding  CUtf's,  TUW's  and  CTUW's,  are  output  in 
side-by-side  decimal  and  octal  format. 
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Program 


INTEGER  UW,SUW,TUW,CTUW,CUW,GUW(73)  ,GCTUW(73)  , 
2MC£g]s3(U6):1,W(   )#rtTEMP(110)  'INCR  (36),  MATCH  (73,  110), 

EQUIVALENCE  (GUW  ( 1)  ,  HCROSS  ( 1)  )  ,  (GCTUW  (1)  ,MCROSS  (7  4)) 
C  INITIALIZATION  AND  PARAMATEE  INPUT 

IDIM=73 

IDIM1=73 

IDIM2=110 
1      READ(5,500)  N,IEPS,IALFA 

HRITE(6#505)  N.IEPS,IALFA 

IE  (N.EC.  (N/2) *2)     GO    TO    5 

K*=  ili<  -\j/2 

LAST=2**K-2 

LASI1=2** (K+1)  -2 

GO    10     10 
5  K=N/2 

LASI=2**  (K+1)-2 

LASI1=2**  (K+2)-2 
10  IZ1-1 

MLAST-N-2 

DO     15    tt=K,MLAST 
15  IZ1=IZH2**M 

IZ2=2**  (N-1)  -2 

MASK=2**N-1 

IFIAG=0 

UW  =  0 

4=1 

1=1 

C  INITIALIZE  EOR  AUTC-C  OF  NEXT  UW 
20     1=0 

SUK=UW 
C  COMMENCE  AUTO-C  LOOP. 

C  THE  FCLIOWING  IN-LINE  SEGMENT  INITIALIZES  AN  N-I  COUNTER, 
C  LOADS  SUW  INTO  A-REGISTER,  SHIFTS  SUW  RIGHT  ONE  BIT, 
C  STORES  RESULTS  IN  ."Sum"  FOR  LATER  USE,  COMPUTES  MOD- 2  SUM 
C  OF  UW  AND  SUW,  COUNTS  NUMBER  OF  ZEROS  IN  FIRST  N-I 
C  POSITIONS  (FROM  RIGHT)  OF  THE  MOD-2  SUM. 

30  1=1+1 
ICCUNT=0 

S  IDA  N 

S  SUE  I 

S  ADE  =077677777 

S  RCH  070100,1 

S  LDA  SUW 

S  SHIFT  C21001 

S  STA  SUW 

S  EOR  UW 

S  SHIFT  0100.1 

S  SKE  =040000000 

S  MPC  ICCUNT 

o        RR  Y  S— 1   1 

C  ARE  WE  INSiSe  APERTURE? 

IF  (I.LE.IALFA)GO  TO  31 
C  ARE  WE  CN  DISTRIBUTION  TAIL? 

IF  (I.GE.N-IEPS)  GO  TO  32 
C  DCES  UW  PASS  AUTO-C  TEST  AT  STEP  I  OUTSIDE  APERTURE? 

IF  (ICCUNT. LT.N-IEPS-I)  GO  TO  30 

GO  TO  35 
C  DOES  UW  PASS  AUTO-C  TEST  AT  STEP  I  INSIDE  APERTURE? 

31  IF  (ICCUNT. LT.N-2*IEPS-I)  GO  TO  30 
GO  TO  35 

DOES  UK  FAIL  ZERO  TEST  ON  DISTRIBUTION  TAIL? 


32     IF  (ICCUNT. NE.O)  GO  TO  35 
C  NOW  ON  DISTRIBUTION  TAIL 
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33 

C 
C 

c 
c 
c 
c 
s 
s 
s 
s 
s 
s 
s 
s 
s 
s 
s 
s 
s 
s 
s 
s 
c 


JFLAG=JFLAG+ 1 

IF  (JFLAG.LT.IEPS)  GO  TO  30 
HAVE  GOOD  UW.   COMPUTE  CUW,  TUW,  CTUW. 

THE  FOLLOWING  IN-LINE  SEGMENT  INITIALIZES  AN  N  COUNTER. 
LOADS  UW  INTO  A-REGISTER,  FORMS  TRANSPOSE  OF  UW  BY  RIGHT 
SHIFTING  ONE  BIT  AT  A  TIME  INTO  B-REGISTER  AND  THEN  CYCLE 
4END-ARCUNL)  SHIFTING  THE  LATTER,  FORMS  COMPLEMENT  OF 
BOTH  UW  AND  TUW  USING  MASK  OF  ONES. 

LDA  N 

ADD  =077677777 

RCH  070100f1 

LDA  UW 

RCH  037710 

SHIFT  C1001 

SHIFT  C16002 

BRX  $-2,\ 

SHIFT  C1 1001 

STB  TUW 

LDA  UW 

EO£  MASK 

STA  CUW 

LDA  TUW 

EOR  MASK 

STA  CTUW 
HAVE  W£  PREVIOUSLY  CHECKED  AND  STORED  THESE 

IF  (IUW. LT.UW.OS.CTUW.LT. UW)  GO  TO  35 
C  DIMENSION  OVERFLOW  CHECK 

IF  (L.GIelDIM)  GO  TO  47 
C  STORE  GCOD  SEQUENCES 

GUW  (L)  =UW 

GCUW  ' 

GT 

GC 

L=L  +  1 
C  PICK  A  NEW  SEQUENCE 
35     IF  (Uw.LT.LA5T)  GO  TO  4  0 

IF  (UW.GE.IZ1)  GO  TO  45 

IF  (UW.LT.LAST1)  GO  TO  52 

J  =  J  +  1 

IF  (N.EQ.  (N/2)*2)  GO  TO  50 

LASI1=2**(K+J) -2 

GO  TO  52 
50     LASTl=2**(K+J+1)-2 
52     LAS1=1AST+2**(K-J) 

UW=UW+3 


SEQUENCES? 


I  W  CI)  =u« 

:uw  (L)=cuw 

:UW(L)=TUW 
:1UH  (L)=CTUW 


55     GO  TO  20 

40     UW=UW+1 

GO  TO  20 
45     IF  (IFLAG.EQ.  1) 

IF1-AG=IFLAG+1 

UW=IZ2 

GO  TO  20 
C  DIMENSION  OVERFLOW 
47     WRITE  (6,506)  UW 

STOP 
C  COMPUTE  NUW  AND  NAUTO 
60     L1=L-1 

NUW=L1 

NAUIO=2*NUW 
C  OUTPUT  GOOD  SEQUENCES 

WRITE  (6,507) 

DO  65  M=1.L" 

WRITE  (6,508 


GO  TO  6  0 


HALT 


(DECIMAL  AND  OCTAL  EORMAT) 


65 

C 

500 
505 
506 

507 

508 


1GTUW(M) ,GCT 


)  M,GUW(M)  .GUW(M) 
UW(M)  ,GCTUW  (M) 


.GCUW(M)  ,GCUW  (M)  ,  GTUW  (M) 


3110) 


FORMAT  ^1«"iN=.',l2,10Xr»IEPS=«,H/10X,'IALFA=',lf) 
FOENAT  «0''« AUTOCORRELATION  SEARCH  TERMINATED  TO  , 
«  P IEVFNT  DIMENSION  OVERFLOW.  LAST  WORD  CHECKED:  •  ,16) 
FOafiTI'i'  • UNIQUE  WORDS  WHICH  PASS  AUTOCORRELATION  , 
t  lEST'//l6x^ia«t25X#»  CUW,  25X, • TU W  ,25X, •  CTUW//) 
FORMAT  <*  •  ,l5,5X,4  (16  ,  3X  ,0  1  0 ,  8X)  ) 
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Cross-CorreJation  Sort 


Variables 


K2   -  Numter  of  sequences  with  which  sequence  beinq 
tested  has  cross-correlated  acceptably.   Used  for 
filling  MATCH  .matrix. 

M2HIGH  -  Highest  number  of  sequences  with  which  any  one 
sequence  has  cross-correlated  acceptably.   Used  tor 
MA1CH  matrix  o-utput.   (This  is  no£  the  maximum  number 
of  cross-correlation  matches  exrstinq.   Correlation  is 
checked,  and  marches  stored  bui;  once,  commencing  with 
the  first  sequence  in  the  AUTO-C  set  and  workinq  up. 
Thus  M2HIGH  does  not  include  matches  between  the 
sequence  beinq  tested  and  previous  sequences  in  the 
AUIO-C  set.) 

Other   variables  are  defined  under  "Arrays"  and 
"Comacn. " 


Output 


(1)  Ihe  UW  and  GTUN  sequences  from  the  AUTO-C  set  which 
are  being  cross-correlated  are  output  in  decimal  format, 
together  with  an  index  number  which  denotes  the  order  in 
whrch  the  sequences  are  tested. 

(2)  The  mat.chinq  matrix  (MATCH)  is  output  in 
two-dimensional  form.  The  elements  of  MATCH  are  the 
index  numbers  corresponding  to  the  AUTO-C  set  of 
sequences.  There  is  one  "row  for  each  AUTO-C  UK  (See 
Section  IV. B).  The  first  element  cf  a  row  is  the  index 
number  of  the  UW,  and  the  remainder  of  the  row  contains 
jhiqher)  index  numbers  of  the  UW's  and  CTUW's  with  which 
it  matches. 
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P  r  pg  r  a  a 


C  FILL    MA1CH  KITH  ZEROS 
DO  103  I=1,IDI.M1 
DO  103  0=1,IDIU2 

103  UA1CH  (IrJ)  =0 
C  PACK  KCROSS 

DO  104  1=1. NUW 
C  OUTPUT  THE  AUTO-C  SET  ELEMENTS  AND  CORRESPONDING  INDICES 

104  MCEG5S  (NUh+I) =GCTUH  (I) 
WRIIE(6,610) 

WRITE  (6,61 1)  (I,MCROSS (I)  ,I=1,NAUIO) 
C  INITIALIZATION 

L=1 

M2  =  2 

M2HIGH=2 
C    SELECT     NEXT    UK    FROil    AUTO-C    SET 

105  UK=MCRCSS 


3FLOH 


DIMENSION  CVERFLOH  TEST 

IF  (L.GT.IDIM1)  GO  TO  155 
STORE  INDEX  NUMBER  OF  UK 


MTEMP  (1)  =L 
INDEX  OF 


C  SET  INDEX  OF  FIRST  AUTO-C  ELEMENT  BEING  CHECKED  AGAINST  UK 

M  =  L  +  1 
C  CHECK  CEOSS-CORRELATION  WITH  ZERO  SHIFT  EETWEEN  UW  AND 
C  SOW.   THE  IN-LINE  SEGMENT  IS  IDENTICAL  IN  PURPOSE  TO  THE 
C  FIRST  IN-IINE  ROUTINE  OF  PREVIOUS  PROGRAM  SEGMENT 
110    ICCUNT-0 

1=0 

JFLAG=0 

S0fc=MCEOSS  (M) 
S      IDA  N 

S-»tvt>  _r\-~i—>r—i-i~i—ii 
HU  U      —  \)   i    I  1}  I    I   I    I   i 

S  RCH  070100,1 

S  LDA  UK 

S  EOR  SUK 

S  SHIEI  01001 

S  SKE  =040000000 

S  MPC  ICCUNT 

C  DCES  PAIR  FAIL  TEST  AT  1=0? 

IF (ICOUNT.GE.N-IEPS)  GO  TO  117 
C  COMMENCE  LCOP  FOR  COMPUTING  RIGHT-SHIFTED  CROSS- 
C  CORRELATION.   THE  IN-LINE  SEGMENT  IS  IDENTICAL  IN  PURPOSE 
C  10  THE  CNE  ABOVE. 
115    1=1+1 

ICOUNT=0 
S      LDA  N 
S      SUE  I 
S      ADC  =077677777 
S      RCH  070100, 1 
S      LDA  SUfc 
S      SHIFT  C21001 
S      STA  SUfc 
S      EOR  UK 
S      SHIFT  01001 
S      SKE  =040000000 
S      MPC  ICCUNT 
c        R  R  X  5—  "^  1 
C  ARE  KE  ON  DISTRIBUTION  TAIL? 

IF  (I.GFoN-IEPS)  GO  TO  140 
C  DOES  PAIR  PASS  TEST  AT  STEP  I? 

IF  (ICOUNT.LT. N-IEPS-I)  GO  TO  115 
C  INCREMENT  INDEX 
117    M=M+1 
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C  HAVE  WE  CHECKED  UW  AGAINST  ALL  HIGHER  ORDER  SEQUENCES? 

IF  (M.LE..NAUTO)  GO  TO  110 
C  FINISHED  WITH  UW.  STORE  MATCHES. 
120    M2=M2-1 

DO  125  K=1,M2 


125   mtemp?k: 


MATCH  (I,K)  =MIEMP(K) 

;)=o 

C  STORE  NUMBER  OF    UW  MATCHES  IN  FIRST  UNUSED  ROW  ELEMENT 

MATCH  (L,IDIM2)  =M2 
C  IS    M2  KIGHFR  THAN  .FOR  ANY  PREVIOUS  UW? 

IF  (M2.GI.M2HIGH)  M2HIGH  =  M2 

M2=2 
C  TIME  TC  GET  NEW  UW,. 

L=L+1 
C  ARE  WE  FINISHED? 

IF  (L.LE.NUW)  GO  TO  105 
C  CROSS-C  SORT  COMPLETE.   OUTPUT  MATCH. 
135    WRITE  (6,655) 

L-L-1 

DC  137  1=1, L 
137    WRIIE(.6,660)  (MATCH  (I ,  J)  ,  J=1  ,  M2HIGH) 

GO  TO  300 
C 

C  GC  TO  MATCHING  ROUTINE 
C 

C  DCES  PAIR  EASS  ZERO  TEST  AT  STEP  I? 
140    IF  (ICOUNT.GT.O)  GO  20  117 

JFLAG=JFLAG+1 
C  ARE  WE  FINISHED  WITH  RIGHT-SHIFI? 

IF  (OFLAG.LT.IEPS)  GO  TO  115 

JFLAG=0 
SUW=MCBOSS (M) 

C  START  LOOP  FOR  COMPUTING  LEFT-SHIFTED  CROSS-CORRELATION. 
C  THE  IN-LINE  SEGMENT  INITIALIZES  N-I_AND_I  ^COUNTERS  ,_  IOADS 
C  SUW  IN'IC  A-iix-GISTE-R/  SHIFTS  5UW  LEFT  ONE  BIT,  STORES  THE 
C  RESULT  FCR  LATER  'USE,  COMPUTES  MOD-2  SUM  OF    UW  AND  SUW, 
C  SHIFTS  THE  MOD-2  SUM  RIGHT  I  BITS  SO  THAT  THE  N  BITS  OF 
C  INTEREST  ARE  RIGHT- JUSTIFIED  IN  THE  A-REGISTER,  COUNTS 
C  NUMBER  OF  ZEROS  IN  FIRST  N-I  POSITIONS  (FROM  RIGHT)  CF 
C  THE  MOE-2  SUM. 
145    1=1+1 

ICCUNT=0 
S      IDA  N 
S      SUE  I 
S      ADD  =077677777 
S      RCH  070100,1 
S      LDA  I 
S      ADD  =077677777 
S      RCH  070100,2 
S      LDA  SUW 
S      SHIFT  C25001 
S      ST A  SUW 
S      ECR  UW 
S      SHIFT  01001 
S      BRX  $-1,2 
S      SHIFT  0100  1 
S      SKE  =040000000 
S      MPC  ICCUNT 
S      BRX  $-3,1 

IF (I.GE.N-IEPS)  GO  TO  150 

IF  (ICOIjNT.LT.  N-IEPS-I)  GO  TO  145 

GO  TO  117  „„_ 

150    IF  (ICOUNT.GT.O)  GO  TO  117 

JFLAG=JFLAG+  I 

IF (JFLAG.LT.IEPS)  GO  TO  145 
C  HAVE  A  GOOD  MATCH.  STORE  IN  MTEMP 

MTEMP  (C2)=M 

M2=M2+1 
C  DIMENSION  OVERFLOW  CHECK 

IF (M2.LT.IDIM2-1)  GO  TO  117 

WRITE  (6,665)  L 
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GO  TO  120 
155    WRITE(6,67Q) 

GO  10  135 
C 

610  FORMAT  ('1'  ,' AUTOCORRELATION  WORDS  FROM  ABOVE  WHICH', 
1'  ARE  EEING  CHECKED  FOR  CROSS-CORRELATION  MATCHES'//) 

611  FORMAT  (6  (1X.I4-*. '  2X,I5,8X) ) 

655    FORMAT  ('0',' GOOD  CROSS-CORRELATION  PAIRS.  THE  NUMBER* 

1'  LISTED  IS  THE  SUBSCRIPT  OF  THE  CORRESPONDING  GCOD', 

2'  AUTOCORRELATION  WORD.'//) 
660    FORMAT  <•<)•  ,  10(1X,  2016/)  ) 
665    FORMAT  ('0' , 'OVERFLOW  OF  DIMENSION  IDIM2  ATTEMPTEL', 

1  '  WHEN  L='  3 1) 
670    FORMAT  ('0' \ 'PROGRAM  TERMINATED  TO  PREVENT  OVERFLOW', 

1»  CF  DIMENSION  IDIM11) 
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Matching  HoMtine 
Var^ ablgs 

HOST  -  length  of  largest  set  found  at  any  given  point  in 
search..   Smaller  sets,  when  found,  are  not  output. 

LTESI  -  Sets  the  last  row  element  which  must  be  looked 
at  based  on  current  value  of  MOST. 

MAT1  -  An  element  from  the  matrix  row  which  has  MAT2  as 
the  leading  element. 

MAT2  -  The  last  element  of  the  current  set  being  formed. 

MAT3  -  An  element  from  the  matrix  row  which  has  one  of 
the  set  members  as  a  leading  element.  Once  it  is 
known  that  a  potential  addition  matches  with  the  first 
and  last  elements  of  the  set,  MAT3  is  used  to  check  if 
it  matches  with  the  other  elements  as  well. 

Other  variables  are  defined  under  "Arrays"  and  "Common." 


Match  sets  are  output  one  line  at  a  time  in  decimal 
format  with  the  autocorrelation  index  preceding  each  element 
Ox.  the  S€  t . 
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IL93E&2 


C  SETS  ABE  FORMED  BY  FINDING  ALL  THOSE  THAT  HAVE  THE  I'TH 

C  BOH  LEADING  ELEMENT  AS  THE  LEADING  ELEMENT  OF  THE  SET. 

C  I 'I  IS  KNOWN  THAT  ALL  ELEMENTS  OF  A  ROW  MATCH  WITH  THE 

C  FIRST  AND  IT  IS  NECESSARY  TO  DETERMINE  WHICH  OF  THESE 

C  IN  TURN  MATCH  WITH  ONE  ANOTHER.   IF  A  ROW  ELEMENT  IS 

C  GREATER  THAN  NUW,  IT  COMES  FROM  THE  CTUW  PORTION  OF  THE 

C  AUTO-C  SET  AND  WILL  NOT  APPEAR  AS  A  ROW  LEADING  ELEMENT. 

C  FOR  THIS  CASE,  THE  CORRESPONDING  UW  ROW  IS  USED  TO 

C  DETERMINE  WHAT  FURTHER  MATCHES  EXIST. 

C 

C 

C  INITIALIZE 

300    WRITE(6,700) 

MOSI=3 

1=1 

M  —  O 

C  START  LOOKING  FOR  SETS  WHICH  HAVE  AUTO-C  NO.  1  AS  LEADING 

C  ELEMENT 

305    K=2 

C  COMPUTE  LIEST 

310    IF  (MATCH  (I, K)  .GT.  NUW)  , LIES T=MATCH  (MATCH  (I,K)  -NUW, 

1IEIK2) +M+1-MOST;  GO  TO  312 

LTEST=MATCH  (MATCH  (I,K)  #IDIM2) + M+ 1-MOST 
C  GO  TO  NEXT  ROW 
312    IF  (MATCH  (I,K)  .GT.O)  GO  TO  315 

1  =  1*1 
C  ARE  WE  FINISHED? 

IF  (I.  GS.  NUW)  STOP 

GO  TO  310 
C  IS  NEXT  ELEMENT  IN  ROW  A  ZERO?   IF  SO,  TIME  TO  GO  TO 
C  NEXT  ROW. 
315    IF  (MATCH  (I,K-r1j  .GT.O)  GO  TO  320 

1  =  1  +  1 

IF  (I.  GT.  NUW)  STOP 

GO  TO  305 
320    L=1 

J  =  K+1 
C  IS    MATRIX  ELEMENT  A    'CTUW'? 
325    IF  (MATCH  (I,K)  .GT.  NUW)  GO  TO  415 
C  START  LOCKING  FOR  .NEXT  ADDITION  TO  MATCH  SET 
327    L=L41 

IF (L.GT.LTESI)  GO  TO  372 

MAI1  =  MATCH  (MATCH  (I,K)  ,  L) 

IF  (MATIoEQoO)  GO    TO  370 
330    MAI2=MATCH  (I,J) 

IF  (MAT2„LEoO(  GO  TO  370 

IF  (MAT2.LE*  MAT/1)  GO  TO  335 

C  DO  WE  HAVE  A  POTENTIAL  ADDITION  TO  SET? 
335    IF  (MATl.EQ„MAT-2)  GO  TO  340 
J  =  J  +  1 

C  IS  THE  CURRENT  SET  LENGTH  EQUAL  TO  TWO?   IF  SO,  TAKE  THE 
C  NEW  ADDITION  WITHOUT  QUESTION. 

C4MUSTICHECKQiF) POTENTIAL  NEW  ADDITION  MATCHES  WITH  OTHER 
C  ELEMENTS  AIREADY  ItN  SET 

M1=M-1 

DO  350  MM=2,M1  n    m^  orn 

IF  (HTEKP(KM)  .GI.NUW)  GO  TO  352 

LL  =  2 
345    MAT3=MATCH  (MTEMP(MM)  ,LL) 
C  IF  FOLLOWING  TtfUE,  NEW  UdITION  NO  GOOD 

IE(MAT3.EQ.,0.OR.MAT3.GT.MAT2)  GO  TO  325 

IF  (MAT3.EQoMAT2)  GO  TO  350 

LL  =  I.L+1 

GO  10  34  5 
350    CONTINUE 
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C  HAVE  NEW  ADDITION  TO  SET 

GO  10  355 
352    MA12=MA12-NUW 

DO  354  MM1  =  MM,J11 

LL  =  2 
353 


GO  TO  325 
.EUoMAT2)  GO  TO  354   " 
LL-IL+1 
GO  10  353 

354  CONTINUE 
MA12=MA12+NUW 

C  STORE  NEW  ADDITION  AND  SET  POINTER 

355  M  =  M  +  1 
MTEMP(M) =MAI2 
INCR  (M)=J 

C  START  LOOKING  EOR  SEQUENCES  WHICH  MATCH  HUH  FIRST  AND 
C  (NEW)  LAST  ELEMENT  OE  SET. 

K=J 

IF  (MATCH  (I,K)  .GT.NUW)  ,LTES1=MATCH  (MATCH  (I,K)  -NUW, 
1IDIK2) +M+1-M0SI;  GO  TO  320 

LTESI=MATCH (MATCH (I,K) ,IDIM2) +M+1-M0ST 

GO  TO  320 
C  STORE  FIRST  TWO  ELEMENTS  OF  2-ELEMENT  SET  AND  SET 
C  POINTERS.   THEN  GO  TO  POINT  WHERE  NEW  ELEMENT  IS  ADD£D 
C  TO  SET 
365    HTEMP(1)=I 

MTEMP  (2j=MATCH,(I,K) 

INCE  (1 

INCfi(2; 

GO  TO 
C 

C  HAVE  EXHAUSTED  POSSIBILITIES  FOR  ADDITIONS  TO  CURREN1  SET 
C 

C  IS  SET  LENGTH  OF  INTEREST?  IF  SO,  OUTPUT;   IN  ANY  CASE, 
C  EACK  Of  TO  NEXT- TO- LAST  ELEMENT  OF  SET  AND  LOOK  FOR 
C  FURTHEE  PC1ENTIAL  ADDITIONS. 
370    If  (M.GE.MOST)  GO  TO  375 
372    IF  (M. EC- 2)  GO  TO  405 

GO  TO  580 
375    MOST=M 

WRITE  (6,705)  (MTEMP (M1)  ,MCROSS  (MTEMP  (M1) )  ,M 1=1, MOST) 
380    K  =  INCR(M-1) 

J^INCR  (M)  +1 

L=1 

IE  (MATCH  (I, K)  .GT.NUW)  GO  TO  450 

MA11=MA1CH  (I,K) 

MA12=MAICH(I,I&CR(M)  ) 
385    L=L+1 

IE  (MATCH  (MAT 1,L)  .LT.MAT2)  GO  TO  335 

H=M-1 

LTEST=MATCH  (MATCH (I,K)  ,IDIM2) +  M+ 1-MOST 

GO  10  327 
405    K=K+1 

IF  (MATCH  (I,K)  .GT.NUW)  ,LTES1=MATCH  (MATCH  (I,K)  -NUW, 
1IDIM2) +M+1-M0ST;  GO  TO  315 

L1ES1=MATCH (MATCH (I,K)  , IDIM2) + M+ 1-MOST 

GO  10  315 
C  REMAINING  STATEMENTS  HANDLE  CHECKS  WHEN  MATRIX  ELEMENT  IN 
C  QUESTION  IS  FROM  CTUW  AUTO-C  SUBSET 
415    MAT3  =  MATCH (I,K) -NUW 
420    L=L+1 

IF  (L.G1.LTEST)  GO  TO  372 

MA11=MA1CH (MAT3,L) 

IF  (MAT1.EQ.0)  GO  TO  370 
425    MA12  =  MA1CH  (I,J)-NUW 

IF  (MAT2.LE.0)  GO  TO  370 

IF  (MAT2„LE.MAT1)  GO  TO  430 

430    IF  (MAT1-EQ.MAT2) , MAT2  =  MAT2  +  NUW ; GO  TO  340 
J  =  J+1 
GO  10  425 
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450 
455 


700 
705 


H(I,K) 

:h(I,in 


-NUW 

CR(M)  )  -NUH 


MAT1=MATCH 

MAT2=MATCi 

1=1  +  1 

IF  (MATCH  (MAT  1,L)  .LT.MAT2)  GO  TO  455 

M=M~1 

LTESI=MATCH  (MA.T  1 ,  IDIM2)  +M+1-MOST 

GO  TO  415 

FORMAT  ('1' ,'GOOE  CROSS-CORRELATICN  SETS'//) 

FORMAT  {'0' ,5  <1X,7  (I4,»)  • , 1 X,I6 , 5X) /)  ) 

END 
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